Ana Luisa Pinheiro 11810407
Ayrton Amaral 11288131
Bruno Groper Morbin 11809875
Caio Febronio 11811482

Instituto de Matemática e Estatística - Universidade de São Paulo | Julho, 2023

Descrição do estudo

Para cada indivíduo pertencente ao estudo, foram amostrados genes individuais e realizada a medição da intensidade da luz em cima de cada um deles. O banco de dados disponibilizará os valores dessa intensidade individualmente para cada gene de cada indivíduo.

1 Problema

Com base no banco de dados, queremos verificar se existe alguma relação entre esses valores da intensidade da luz em cima de cada gene com características pessoais de cada indivíduo.

2 Proposta

Para a resolução do problema, agruparemos os genes pelo valor da intensidade de luz e utilizaremos a informação mútua vista em aula para checarmos se os genes introduzidos em cada cluster possuem distribuição parecida ou não. Após essa checagem, será feita uma distribuição única da intensidade de luz para cada cluster e a mesma será utilizada para caracterizar cada conjunto de cluster em relação às características do indivíduo.

2.1 Roteiro

  • Análise descritiva de todos indivíduos do estudo;
  • Sintetização do conjunto de genes que apresentam intensidade de luz semelhantes pareadas por indivíduos, agrupando por similaridade através de técnica de agrupamento hierárquico;
  • Observar a consistência entre os genes agrupados e as intensidades de luz apresentada para cada indivíduo;
  • Unificar a distribuição de intensidade de luz contida em cada grupo de genes;
  • Identificar se há possibilidade de discriminar os níveis de cada variável característica do conjunto de indivíduos, podendo concluir quais grupos de genes têm relação com cada característica.

2.2 Código

# Carregando pacotes
library(tidyverse)
library(dplyr)
library(cluster)
library(infotheo)
library(stringr)
library(patchwork)

2.2.1 Conjunto de dados

load("glioma.RData") # geneInfo ; gliomaGSE52009 ; targetInfoGlioma
glioma <- gliomaGSE52009; as.data.frame(glioma)
info <- targetInfoGlioma; rownames(info) <- NULL; info |> select(colnames(info[,-1]),FileName)
str(info[,-1]) # ignorando a coluna FileName
'data.frame':   120 obs. of  7 variables:
 $ gender      : chr  "female" "female" "female" "unknown" ...
 $ age         : num  37 61 37 -100 31 42 42 38 32 40 ...
 $ geoAccession: chr  "GSM1257398" "GSM1257399" "GSM1257400" "GSM1257401" ...
 $ sampleInfo  : chr  "Astrocytoma.a-0253" "Astrocytoma.a-0258" "Astrocytoma.a-0281" "Astrocytoma.a-0285" ...
 $ diagnostic  : chr  "astrocytoma" "astrocytoma" "astrocytoma" "astrocytoma" ...
 $ datasetId   : chr  "Glioma52009" "Glioma52009" "Glioma52009" "Glioma52009" ...
 $ tissue      : chr  "brain" "brain" "brain" "brain" ...
info$gender <- factor(info$gender); levels(info$gender)
 "female"  "male"    "unknown"
ggplot(info, aes(x=reorder(gender, -table(gender)[gender])))+
  geom_bar(aes(fill=gender), color="transparent")+
  scale_fill_manual(values=c(male="#3E67A3",female="#A34336",unknown="#7F8F85"))+
  geom_text(stat = 'count', aes(label = paste0(round((after_stat(count)/sum(after_stat(count)))*100), "%")), vjust = -0.5) +
  scale_y_continuous(expand = expansion(mult=c(0,.2)))+
  labs(x=NULL,y=NULL, title="Gênero")+
  guides(fill="none")

info$diagnostic <- factor(info$diagnostic); levels(info$diagnostic)
 "anaplastic.astrocytoma"       "anaplastic.oligodendrocytoma"
 "anaplastic.oligodendroglioma" "astrocytoma"                 
 "glioblastoma"                 "oligodendroglioma"           
ggplot(info, aes(x=reorder(diagnostic, -table(diagnostic)[diagnostic])))+
  geom_bar(aes(fill=diagnostic), color="transparent")+
  scale_fill_grey(end = 0.9, start=.5)+
  geom_text(stat = 'count', aes(label = paste0(round((after_stat(count)/sum(after_stat(count)))*100), "%")), vjust = -0.5) +
  scale_y_continuous(expand = expansion(mult=c(0,.15)))+
  labs(x=NULL,y=NULL, title="Diagnóstico")+
  theme(axis.text.x = element_text(angle = 45,hjust = 1, vjust=1))+
  guides(fill="none")

range(info$age)
 -100   70
cat(paste0(sum(info$age<=0)," entradas inválidas para idade  => ",round(sum(info$age<=0)/nrow(info)*100,2), "% da amostra")) # Porcentagem de dados inválidos para idade
17 entradas inválidas para idade  => 14.17% da amostra
range(info[which(info$age>0),]$age)
 17 70
ggplot(subset(info[which(info$age>0),]), aes(x = age, y = after_stat(density))) +
  geom_histogram(aes(y = ..density..), fill = "skyblue", color = "#0c0c0c", binwidth = 5, alpha = 0.9) +
  geom_density(color = "cyan", linetype = "solid", linewidth = 1, fill="transparent") +
  labs(x = "Idade", y = "Frequência relativa", subtitle ="(somente dados válidos)", title="Distribuição de idade na amostra")+
  scale_y_continuous(expand = expansion(mult=c(0,.20)))+
  scale_x_continuous(n.breaks = 20)

unique(info$datasetId)
 "Glioma52009"
unique(info$tissue)
 "brain"

2.2.2 Agrupamento dos genes

Aplica-se um agrupamento hierárquico para simplificar conjunto de genes que apresentam mesma intensidade de luminosidade para cada indivíduo. Dessa forma, e levando em consideração que os dados são coletados sob a mesma medida, opta-se por não padronizá-los. Caso fossem padronizados, algumas intensidades não expressivas poderiam ser conectadas com outras mais expressivas originalmente (quando considerados os valores não padronizados).

d <- dist(glioma, method = "euclidean")

Para a junção dos clusters, tem-se algumas ordens de hierarquização, tais quais:

Linkages representação

  • Single linkage: É a distância mais curta entre quaisquer dois pontos nos dois grupos.

  • Complete linkage: É o oposto do single linkage. É a distância mais longa entre quaisquer dois pontos nos dois grupos.

  • Average linkage: É a média das distâncias entre cada ponto de um grupo para todos os pontos do outro grupo.

  • Centroid linkage: A distância entre o ponto central de um grupo ao ponto central do outro grupo.

  • Ward’s linkage: Uma combinação dos métodos average e centroid. A variância dentro do grupo é calculada determinando o ponto central do grupo e a distância das observações em relação ao centro. Ao tentar mesclar dois grupos, a variância é calculada entre os grupos e os grupos são mesclados se sua variância for menor em comparação com outras combinações.

hc<-hclust(d, method = "ward.D") # Hierarchical Cluster com linkage por Ward.D
num_cluster_init <- 20; init_clusters = cutree(hc, num_cluster_init) # pegando os clusters iniciais 
# Exemplo do um cluster selecionado inicialmente
as.data.frame(glioma[which(init_clusters==1),])
genes.clusters <- list()
for(clust in 1:num_cluster_init){
  gp <- glioma[which(init_clusters==clust),] # selecionando o cluster
  
  if(is.null(nrow(gp))) {
    next
  } else{
    bounds <- apply(gp,MARGIN = 2, FUN = function(x) quantile(x,c(0.15,0.85))) # definindo os intervalos interquantil para cada indivíduo nesse conjunto de genes
    
    # identificando cada intensidade dentro do cluster se está dentro do interquantil do indíviduo
    in_bound <- as.data.frame(lapply(1:ncol(gp), FUN= function(i) {as.numeric(between(gp[,i],bounds[1,i],bounds[2,i]))})); colnames(in_bound) <- colnames(bounds); rownames(in_bound) <- rownames(gp)
    
    # manter no cluster apenas os genes que se apresentaram bastante no todo da amostra dentro do intenquartil da intensidade de cada indivíduo 
    select_genes <- apply(in_bound, MARGIN = 1, function(x) sum(x)) >= 0.7*ncol(glioma)
    select_names <- names(which(select_genes==T))
    if(length(select_names)>1){
      genes.clusters <- append(genes.clusters,list(select_names))
    }
  }
}

for(i in 1:length(genes.clusters)){
  if(i==1) cat("--- Genes:\n")
  cat(paste0("\nCluster ",i,":\n"))
  cat(genes.clusters[[i]],sep=" | ")
  cat("\n")
}
--- Genes:

Cluster 1:
1 | 1001 | 100128288 | 10036 | 10042 | 100505741 | 10096 | 10188 | 10238 | 10301 | 10312 | 10362 | 10391 | 10434 | 10459 | 1046 | 10466 | 10590 | 10595 | 10603 | 10865 | 10885 | 10970 | 11000 | 11004 | 11014 | 11068 | 11117 | 11237 | 112616 | 11277 | 11285 | 11318 | 11320 | 1138 | 115201 | 117143 | 126133 | 126789 | 127002 | 128637 | 129138 | 129303 | 144233 | 144568 | 145389 | 1455 | 146562 | 146853 | 148709 | 1513 | 152137 | 1535 | 157247 | 1586 | 1601 | 162073 | 162417 | 162515 | 1643 | 165140 | 166012 | 1737 | 1846 | 1912 | 1992 | 2030 | 2034 | 204851 | 205564 | 2064 | 2180 | 219293 | 219855 | 221035 | 2235 | 22838 | 22916 | 2296 | 2303 | 2305 | 2307 | 23208 | 2326 | 23304 | 23597 | 23659 | 23765 | 245812 | 254359 | 254552 | 255031 | 25825 | 25853 | 25865 | 259240 | 26037 | 26531 | 26578 | 2677 | 27143 | 27349 | 2802 | 2815 | 2829 | 282997 | 283219 | 283298 | 283933 | 284119 | 284266 | 284346 | 285877 | 286053 | 286102 | 286133 | 286826 | 2875 | 2882 | 28951 | 28959 | 29108 | 29886 | 29928 | 2997 | 307 | 30844 | 3216 | 3219 | 3273 | 3385 | 339122 | 340351 | 3430 | 3460 | 349149 | 350383 | 3695 | 3703 | 374907 | 374969 | 375775 | 388759 | 388815 | 390883 | 3959 | 4017 | 4047 | 4054 | 410 | 4161 | 4210 | 4258 | 4318 | 4598 | 4627 | 47 | 474344 | 4839 | 4864 | 4927 | 4952 | 5002 | 5020 | 50848 | 51225 | 51388 | 51393 | 5154 | 5175 | 51773 | 5211 | 5266 | 5292 | 5322 | 5359 | 5424 | 5430 | 54521 | 54538 | 54555 | 54815 | 54856 | 54919 | 54962 | 54964 | 54967 | 54995 | 55039 | 55127 | 55178 | 55274 | 55276 | 55287 | 55471 | 55536 | 55599 | 55692 | 55707 | 55775 | 55816 | 55858 | 5598 | 5606 | 56122 | 56245 | 5630 | 56731 | 56981 | 56983 | 57124 | 57125 | 57190 | 57232 | 57690 | 5783 | 5819 | 58485 | 5898 | 58986 | 5921 | 59274 | 59342 | 598 | 60559 | 6184 | 6195 | 6282 | 6329 | 63901 | 64215 | 64427 | 64599 | 6464 | 6503 | 6574 | 65996 | 6628 | 6645 | 6753 | 6915 | 696 | 705 | 7059 | 7096 | 7106 | 7132 | 7252 | 7326 | 750 | 775 | 7791 | 7802 | 79080 | 79148 | 79324 | 79567 | 79686 | 79701 | 79759 | 79784 | 79845 | 79901 | 80004 | 80199 | 80271 | 80741 | 80761 | 80776 | 81872 | 83596 | 83743 | 84279 | 84318 | 8473 | 84941 | 84948 | 85007 | 85458 | 857 | 8576 | 8764 | 8773 | 8908 | 89848 | 90187 | 90417 | 9047 | 9064 | 90861 | 91012 | 9136 | 91574 | 91748 | 92241 | 93145 | 94056 | 9473 | 9702 | 9749 | 976 | 9918 | 9938

Cluster 2:
10001 | 100130815 | 100132319 | 10076 | 10107 | 10219 | 10246 | 10316 | 10336 | 10380 | 10564 | 10568 | 1071 | 10827 | 11128 | 11149 | 11184 | 11211 | 11227 | 11238 | 112479 | 112574 | 114625 | 114825 | 115004 | 1161 | 116519 | 118460 | 1187 | 118788 | 119 | 120379 | 120534 | 1215 | 122616 | 122618 | 123263 | 1244 | 124599 | 124930 | 124975 | 125919 | 125958 | 126123 | 126626 | 127391 | 130589 | 132141 | 132671 | 133308 | 134265 | 134285 | 134430 | 1360 | 1365 | 138429 | 1385 | 139285 | 139599 | 139818 | 1417 | 142683 | 143098 | 143425 | 143503 | 144097 | 144100 | 144321 | 144811 | 145858 | 146279 | 146456 | 146850 | 147372 | 147429 | 147719 | 147837 | 148109 | 148229 | 1489 | 149830 | 150248 | 1511 | 151742 | 152815 | 153478 | 153642 | 1572 | 157657 | 157773 | 158747 | 158787 | 1604 | 161779 | 163049 | 163115 | 163590 | 163747 | 163859 | 163882 | 1644 | 165082 | 166378 | 167127 | 170393 | 170679 | 1733 | 1820 | 1907 | 192134 | 1943 | 1990 | 200014 | 200197 | 200205 | 2017 | 202309 | 204962 | 2103 | 2113 | 2124 | 2157 | 219743 | 220032 | 220979 | 221092 | 221785 | 222662 | 2251 | 2260 | 2280 | 22801 | 22888 | 2299 | 23023 | 2304 | 23078 | 23090 | 23094 | 23137 | 23143 | 23144 | 23195 | 2321 | 2322 | 23225 | 23266 | 23276 | 23306 | 23409 | 23476 | 23533 | 23548 | 23626 | 23650 | 23742 | 23786 | 239 | 24146 | 2444 | 245 | 246269 | 246329 | 248 | 2538 | 253943 | 255022 | 255725 | 255738 | 256643 | 2568 | 257106 | 257160 | 257629 | 258 | 25820 | 25895 | 25896 | 25914 | 26007 | 26036 | 26083 | 26093 | 26136 | 26184 | 26205 | 26235 | 26273 | 26280 | 26287 | 2673 | 26735 | 27005 | 27039 | 27042 | 27091 | 27097 | 2710 | 27125 | 27235 | 27283 | 27297 | 2760 | 2797 | 28234 | 2827 | 2830 | 283130 | 283464 | 283600 | 283777 | 284 | 284076 | 284086 | 284098 | 284099 | 284110 | 285525 | 285659 | 286234 | 286464 | 2866 | 2877 | 2898 | 2948 | 2957 | 29767 | 29850 | 29904 | 29907 | 29966 | 29986 | 30009 | 3044 | 3045 | 312 | 3175 | 3196 | 3239 | 326 | 3305 | 338323 | 339221 | 339240 | 339488 | 339906 | 3425 | 343 | 344658 | 344752 | 347365 | 3554 | 3575 | 359821 | 3660 | 3683 | 374864 | 387263 | 388323 | 388581 | 392138 | 3976 | 3980 | 3993 | 399979 | 400673 | 401232 | 403315 | 4041 | 405 | 4113 | 4218 | 4238 | 4277 | 4286 | 4300 | 4317 | 4337 | 4361 | 440503 | 442194 | 4440 | 4524 | 4589 | 4647 | 4668 | 4680 | 4863 | 487 | 4881 | 4907 | 4939 | 4988 | 5004 | 5049 | 50618 | 5069 | 50846 | 5089 | 51091 | 51174 | 51200 | 51207 | 51276 | 51362 | 51378 | 51400 | 51439 | 51449 | 51471 | 51545 | 51666 | 51702 | 51765 | 51804 | 5304 | 5366 | 53820 | 53828 | 53829 | 53981 | 54545 | 54581 | 54583 | 54625 | 5465 | 54660 | 54718 | 54726 | 54738 | 54758 | 54777 | 54831 | 54861 | 54947 | 54979 | 55024 | 55051 | 55102 | 55132 | 55142 | 55164 | 55184 | 55223 | 55243 | 55266 | 55285 | 55334 | 5534 | 55435 | 55620 | 55628 | 55662 | 55728 | 55761 | 55814 | 55869 | 55879 | 5599 | 56172 | 5623 | 56269 | 56666 | 56834 | 5697 | 57010 | 57085 | 57115 | 57143 | 57169 | 5733 | 57465 | 57504 | 57538 | 57621 | 57674 | 57679 | 5780 | 57826 | 5810 | 5820 | 58495 | 5873 | 5913 | 5916 | 5922 | 60343 | 6092 | 6097 | 6123 | 6199 | 6296 | 6314 | 63910 | 63948 | 63951 | 64081 | 64150 | 64282 | 64288 | 64375 | 645644 | 645974 | 64641 | 646513 | 64689 | 64708 | 64783 | 64841 | 649 | 6490 | 64927 | 65258 | 65267 | 6527 | 652995 | 6534 | 6561 | 6580 | 65988 | 6646 | 6667 | 668 | 6773 | 6776 | 6789 | 6847 | 6892 | 6911 | 6928 | 7007 | 7032 | 7036 | 7104 | 7112 | 7145 | 7172 | 7182 | 7185 | 7221 | 7291 | 729975 | 7358 | 7369 | 7373 | 740 | 7469 | 7481 | 7547 | 7554 | 7565 | 7629 | 7633 | 766 | 7705 | 7712 | 7780 | 7813 | 79018 | 79066 | 79368 | 79411 | 79673 | 79836 | 79865 | 79879 | 79959 | 79987 | 8001 | 80014 | 80055 | 80059 | 80071 | 80119 | 80152 | 80153 | 80177 | 80216 | 80307 | 80709 | 80765 | 80824 | 80854 | 8139 | 8161 | 81627 | 81629 | 81789 | 81796 | 81832 | 81857 | 8313 | 83463 | 83639 | 83694 | 83695 | 83734 | 83795 | 8396 | 83998 | 84033 | 84076 | 84107 | 84135 | 84330 | 8440 | 84536 | 84548 | 84612 | 8462 | 84787 | 84793 | 84836 | 84874 | 84895 | 8501 | 8521 | 85291 | 85363 | 8540 | 85456 | 8547 | 8558 | 8658 | 8676 | 8800 | 8840 | 8846 | 8867 | 8894 | 8915 | 89782 | 89910 | 89978 | 90025 | 90102 | 90246 | 90407 | 9063 | 90693 | 90957 | 91351 | 9139 | 91544 | 9159 | 91734 | 9187 | 92014 | 92104 | 93134 | 9342 | 93649 | 9370 | 94086 | 9420 | 9436 | 9439 | 944 | 9468 | 9498 | 9537 | 9542 | 9623 | 9649 | 9655 | 96764 | 9689 | 9696 | 9818 | 9841 | 9855 | 9910 | 9982

Cluster 3:
10011 | 100130175 | 100130856 | 10015 | 10020 | 10048 | 10067 | 10075 | 10078 | 10081 | 10101 | 10103 | 10147 | 10152 | 10153 | 10195 | 10237 | 10277 | 10291 | 103 | 10327 | 10445 | 10469 | 10478 | 10482 | 10483 | 10557 | 10591 | 10616 | 10620 | 10621 | 10634 | 10638 | 10652 | 10659 | 10667 | 10712 | 10725 | 10730 | 10743 | 10772 | 10845 | 10914 | 10945 | 11047 | 11051 | 11066 | 11072 | 11104 | 11112 | 11161 | 11272 | 112936 | 113829 | 114926 | 115817 | 116115 | 116461 | 1180 | 118472 | 1185 | 119559 | 120227 | 121642 | 122704 | 123879 | 124152 | 124491 | 124565 | 124936 | 124976 | 127018 | 127829 | 128 | 1291 | 129531 | 129831 | 130617 | 132 | 132014 | 133957 | 134429 | 135138 | 138065 | 1395 | 140691 | 142678 | 1432 | 144348 | 1452 | 1457 | 146223 | 147179 | 147650 | 147807 | 1479 | 147906 | 148022 | 149076 | 149951 | 150383 | 150678 | 150684 | 151516 | 153396 | 158055 | 159195 | 160298 | 1616 | 163 | 168544 | 170622 | 170680 | 1760 | 1762 | 1785 | 1802 | 192286 | 1949 | 196403 | 197258 | 201627 | 203 | 203245 | 2071 | 2110 | 2117 | 2186 | 219793 | 220002 | 222068 | 222256 | 22796 | 22803 | 22818 | 22845 | 22846 | 22879 | 22880 | 22882 | 22897 | 22919 | 23042 | 23060 | 23064 | 23082 | 23155 | 2317 | 23185 | 23193 | 23198 | 23248 | 23288 | 23291 | 23317 | 23334 | 23344 | 23367 | 23382 | 23387 | 23401 | 23435 | 23469 | 23506 | 23515 | 23523 | 23524 | 23595 | 23608 | 23613 | 23635 | 23647 | 23788 | 2561 | 257068 | 257236 | 25771 | 25792 | 25799 | 2580 | 25807 | 2585 | 25855 | 25870 | 25904 | 25911 | 25950 | 25979 | 25980 | 25988 | 26001 | 26088 | 26156 | 261726 | 26260 | 26268 | 26269 | 26505 | 26528 | 27013 | 27301 | 27434 | 283237 | 283899 | 284827 | 285971 | 2870 | 29035 | 29115 | 29126 | 2959 | 29855 | 29883 | 29954 | 29978 | 3064 | 30968 | 3192 | 323 | 329 | 337867 | 345651 | 348180 | 353116 | 3588 | 3622 | 3720 | 3843 | 3850 | 3880 | 3882 | 388272 | 392 | 400569 | 4059 | 4140 | 4204 | 4215 | 4247 | 4331 | 4436 | 4482 | 4580 | 4790 | 4801 | 4817 | 493856 | 4968 | 5018 | 5048 | 50650 | 50813 | 5082 | 50835 | 51076 | 51093 | 51118 | 51182 | 51250 | 51251 | 51271 | 5134 | 51397 | 5147 | 51526 | 51552 | 51594 | 51603 | 51619 | 51649 | 51663 | 5170 | 51728 | 51734 | 51754 | 51762 | 51763 | 5191 | 5252 | 5253 | 5310 | 5356 | 540 | 54148 | 5422 | 5428 | 54431 | 54432 | 54455 | 54471 | 54494 | 54780 | 54808 | 54812 | 54840 | 54868 | 54878 | 54994 | 55000 | 55003 | 55020 | 55023 | 55035 | 55069 | 55070 | 55072 | 55114 | 55122 | 55131 | 55145 | 55152 | 55156 | 5520 | 55226 | 55245 | 55251 | 55260 | 5533 | 55374 | 5550 | 55593 | 55596 | 55622 | 55654 | 55663 | 55667 | 55683 | 55687 | 55723 | 55731 | 55752 | 55755 | 55759 | 55763 | 5580 | 55829 | 55844 | 55850 | 55904 | 55957 | 56061 | 56181 | 56259 | 56897 | 56912 | 56916 | 56926 | 56941 | 56984 | 56994 | 56996 | 57001 | 57019 | 57140 | 57148 | 57149 | 57187 | 57488 | 57505 | 57551 | 57610 | 57643 | 57645 | 57658 | 57665 | 5822 | 5826 | 58486 | 585 | 5887 | 5889 | 5927 | 5970 | 5981 | 6039 | 60506 | 60673 | 60686 | 6256 | 63027 | 6342 | 63906 | 63917 | 6392 | 64083 | 64118 | 64175 | 64218 | 642852 | 6432 | 64756 | 64760 | 64772 | 64786 | 64969 | 64979 | 65057 | 6597 | 65982 | 6605 | 6626 | 670 | 6714 | 6728 | 6730 | 6731 | 6734 | 6811 | 6836 | 6868 | 6882 | 6885 | 6903 | 6949 | 7016 | 7158 | 717 | 7174 | 7249 | 7265 | 729614 | 7322 | 7327 | 7337 | 7375 | 7514 | 7536 | 7756 | 7841 | 7866 | 7879 | 78991 | 79029 | 79034 | 79078 | 79089 | 79090 | 79183 | 79188 | 79568 | 79585 | 79596 | 79603 | 79624 | 79683 | 79712 | 79714 | 79726 | 79847 | 79871 | 7988 | 79891 | 79893 | 79894 | 80010 | 80011 | 80025 | 80218 | 80232 | 80267 | 80305 | 80315 | 80335 | 80727 | 80737 | 80746 | 80762 | 8086 | 81034 | 813 | 81309 | 81328 | 81558 | 81559 | 81562 | 81566 | 81858 | 819 | 8195 | 8204 | 832 | 8334 | 83480 | 83590 | 83638 | 83692 | 83696 | 83759 | 83862 | 8398 | 84062 | 84065 | 84124 | 84132 | 84138 | 84167 | 84225 | 8428 | 84284 | 84294 | 84311 | 84336 | 84340 | 84518 | 847 | 84706 | 84792 | 84816 | 84866 | 84916 | 84925 | 84942 | 84946 | 84950 | 84954 | 8498 | 85025 | 8517 | 85441 | 8545 | 8559 | 8560 | 8562 | 8568 | 8570 | 8704 | 8723 | 8798 | 8818 | 8824 | 8893 | 8897 | 8924 | 89766 | 90024 | 9013 | 90580 | 90806 | 9091 | 9122 | 9128 | 91300 | 91612 | 91647 | 9166 | 91663 | 91695 | 91746 | 9188 | 92105 | 92154 | 9217 | 9246 | 92597 | 9278 | 9315 | 93190 | 93210 | 9330 | 9337 | 9392 | 9406 | 94101 | 94103 | 9453 | 9491 | 95 | 9519 | 9521 | 954 | 9541 | 9552 | 9557 | 9643 | 9656 | 9657 | 973 | 9732 | 9761 | 9777 | 9784 | 9794 | 9797 | 9862 | 9885 | 9894 | 9903 | 9907 | 9917 | 9919 | 9925 | 9933

Cluster 4:
100129503 | 10017 | 10019 | 100192379 | 10022 | 1003 | 1004 | 10050 | 10082 | 1010 | 10125 | 10149 | 10168 | 10198 | 10242 | 10272 | 10311 | 10390 | 10406 | 10462 | 1048 | 10556 | 10562 | 10584 | 10587 | 10607 | 10610 | 10640 | 10665 | 10677 | 10721 | 10747 | 10750 | 10753 | 10786 | 10800 | 1081 | 10847 | 10861 | 10864 | 10867 | 1087 | 10953 | 11009 | 1103 | 11043 | 11045 | 1107 | 11136 | 11147 | 11181 | 1121 | 11213 | 112401 | 11245 | 11276 | 11279 | 11317 | 11330 | 11338 | 113510 | 1142 | 114609 | 114771 | 114897 | 115560 | 115572 | 115795 | 115992 | 116179 | 116349 | 116969 | 120935 | 121457 | 122769 | 122809 | 122830 | 122876 | 123036 | 124404 | 124535 | 124583 | 124857 | 125488 | 126006 | 126014 | 126375 | 1269 | 127544 | 127670 | 127833 | 127845 | 128486 | 128497 | 129450 | 129530 | 129563 | 1303 | 1305 | 1306 | 130752 | 131566 | 131873 | 133746 | 134728 | 136371 | 137835 | 138050 | 1386 | 139324 | 139716 | 139886 | 140460 | 140609 | 140699 | 140738 | 140850 | 140881 | 140947 | 1414 | 1431 | 143379 | 1435 | 143501 | 143686 | 1440 | 1441 | 145508 | 147081 | 147798 | 147945 | 148137 | 148789 | 149233 | 149483 | 149708 | 149840 | 1500 | 150000 | 150771 | 151056 | 151306 | 151636 | 152816 | 152831 | 152992 | 153201 | 153241 | 153769 | 154089 | 1543 | 155368 | 157769 | 158506 | 1588 | 159090 | 1591 | 1609 | 162333 | 162394 | 162655 | 163154 | 164153 | 165545 | 166785 | 167465 | 169522 | 169693 | 170384 | 171177 | 1776 | 1794 | 1800 | 1806 | 1807 | 181 | 1812 | 1814 | 1815 | 1816 | 1836 | 1889 | 1894 | 1896 | 1903 | 192666 | 1953 | 1962 | 196951 | 197259 | 1991 | 199720 | 199786 | 199964 | 199974 | 200081 | 200909 | 201176 | 201232 | 201283 | 201294 | 201625 | 201626 | 202051 | 202181 | 203430 | 2078 | 2118 | 2121 | 2158 | 219473 | 219995 | 2204 | 220972 | 220992 | 221264 | 2213 | 221472 | 221656 | 222255 | 222484 | 222545 | 2244 | 2249 | 225689 | 2277 | 22802 | 22847 | 22854 | 22900 | 22915 | 22943 | 22999 | 23007 | 23080 | 23085 | 23093 | 23162 | 23181 | 23224 | 23235 | 2324 | 23247 | 23252 | 23275 | 2328 | 2330 | 2334 | 23498 | 23529 | 23567 | 23581 | 23607 | 23682 | 24145 | 245936 | 246175 | 249 | 2526 | 2529 | 252983 | 253012 | 253128 | 2532 | 253827 | 2539 | 253980 | 254013 | 254042 | 254268 | 254427 | 255027 | 2551 | 255189 | 2553 | 255877 | 256076 | 256144 | 256302 | 256646 | 257 | 257101 | 257397 | 25803 | 25819 | 25833 | 25884 | 25885 | 25888 | 25928 | 25937 | 25961 | 25974 | 25981 | 26 | 260436 | 26190 | 2620 | 2622 | 26254 | 26261 | 26262 | 26301 | 26576 | 2658 | 26659 | 266747 | 26747 | 2694 | 27 | 27019 | 27036 | 27071 | 2712 | 27123 | 27124 | 27127 | 27178 | 27180 | 2719 | 27240 | 27291 | 27293 | 27304 | 27314 | 27334 | 2740 | 2762 | 2784 | 2793 | 2799 | 2813 | 2826 | 282618 | 283078 | 283197 | 2832 | 283316 | 283349 | 283588 | 283849 | 283932 | 283948 | 284114 | 284185 | 284242 | 284349 | 285051 | 285282 | 2853 | 285512 | 285668 | 285672 | 286006 | 286410 | 286527 | 286753 | 28964 | 29063 | 29075 | 2909 | 2912 | 2967 | 29777 | 2978 | 2980 | 29800 | 29842 | 29851 | 29881 | 29924 | 29926 | 29930 | 29933 | 2994 | 29948 | 29949 | 30000 | 3003 | 3037 | 3046 | 3050 | 3054 | 3062 | 3067 | 3070 | 3081 | 3101 | 317671 | 320 | 3214 | 3266 | 3284 | 333929 | 3352 | 3357 | 337879 | 337959 | 3386 | 338651 | 339059 | 339451 | 339456 | 339855 | 340168 | 340205 | 340393 | 341276 | 344758 | 346606 | 3480 | 348093 | 348825 | 352909 | 352999 | 353219 | 3557 | 3568 | 3589 | 3604 | 3656 | 3664 | 3674 | 3690 | 3712 | 3742 | 374354 | 375035 | 375287 | 375341 | 3758 | 3772 | 377677 | 3777 | 3818 | 3820 | 3832 | 3836 | 3860 | 387129 | 387496 | 387694 | 387700 | 388121 | 388882 | 389123 | 389136 | 389690 | 390181 | 3957 | 3973 | 399 | 3998 | 401010 | 401105 | 401236 | 4015 | 4063 | 4070 | 4091 | 4117 | 412 | 4143 | 4160 | 417 | 4194 | 420 | 4283 | 4299 | 4311 | 4321 | 440068 | 440072 | 440101 | 440738 | 440955 | 441295 | 4486 | 4606 | 4628 | 4641 | 4671 | 474354 | 4771 | 4798 | 4818 | 4855 | 4862 | 4899 | 4917 | 4919 | 4920 | 4982 | 4986 | 5023 | 5071 | 5074 | 5081 | 5087 | 51005 | 51098 | 51127 | 51148 | 51194 | 5130 | 51312 | 51332 | 51350 | 5140 | 51409 | 51438 | 5144 | 51535 | 51550 | 51561 | 51585 | 5165 | 51684 | 51744 | 51778 | 518 | 5284 | 5293 | 5296 | 5308 | 53354 | 5338 | 53405 | 53615 | 53833 | 54097 | 54101 | 5444 | 54457 | 54474 | 54502 | 5451 | 5453 | 54677 | 54796 | 54798 | 54807 | 54814 | 54822 | 54834 | 54845 | 54853 | 54869 | 54879 | 54880 | 54905 | 54993 | 55025 | 55105 | 55106 | 55124 | 55135 | 55163 | 5519 | 55193 | 552 | 55236 | 55270 | 55293 | 55308 | 55320 | 554 | 554202 | 55432 | 55472 | 55500 | 55507 | 55558 | 55567 | 55568 | 55586 | 55589 | 55632 | 55655 | 55721 | 5578 | 55783 | 5582 | 5603 | 5607 | 56169 | 56244 | 56253 | 5629 | 56311 | 56474 | 56479 | 566 | 56606 | 56667 | 56886 | 56893 | 56923 | 56929 | 56940 | 56949 | 56957 | 57053 | 57089 | 57116 | 57119 | 57154 | 57182 | 57228 | 5731 | 57393 | 5740 | 57402 | 57404 | 57459 | 57462 | 57464 | 57529 | 5753 | 57539 | 57589 | 57650 | 57727 | 57818 | 57823 | 5795 | 5802 | 58496 | 58498 | 58503 | 58510 | 58513 | 586 | 5881 | 5918 | 59271 | 59284 | 5966 | 5991 | 6010 | 60437 | 60529 | 60558 | 607 | 623 | 6343 | 6361 | 6374 | 63971 | 63977 | 64005 | 64091 | 64122 | 64127 | 64174 | 64220 | 6444 | 644809 | 645431 | 64579 | 6462 | 646658 | 64802 | 64848 | 65243 | 65244 | 65266 | 6532 | 6559 | 6570 | 658 | 66037 | 6654 | 6699 | 6710 | 6747 | 6777 | 683 | 6839 | 686 | 6877 | 688 | 6886 | 6932 | 6934 | 7022 | 7052 | 7054 | 7066 | 7124 | 7128 | 7142 | 7157 | 722 | 726 | 728175 | 728743 | 7289 | 7292 | 729420 | 730101 | 7380 | 7428 | 7579 | 7584 | 7596 | 7634 | 7678 | 7739 | 7753 | 7761 | 7767 | 778 | 7784 | 781 | 7903 | 79075 | 79098 | 79170 | 79369 | 796 | 79611 | 79646 | 79648 | 79696 | 79742 | 79788 | 79802 | 79807 | 79814 | 79839 | 79857 | 79861 | 79866 | 79872 | 79917 | 79930 | 79934 | 79946 | 79947 | 79974 | 79979 | 80069 | 80110 | 80128 | 80143 | 80167 | 80176 | 80204 | 80223 | 80237 | 8029 | 80326 | 8038 | 80712 | 80817 | 80821 | 81492 | 81605 | 81691 | 81786 | 81833 | 81839 | 8292 | 8315 | 83401 | 83416 | 83444 | 83445 | 83468 | 83483 | 83544 | 8359 | 83714 | 83729 | 8383 | 83854 | 83856 | 84068 | 84069 | 84077 | 84102 | 8416 | 84182 | 84188 | 84190 | 84197 | 84243 | 84267 | 84376 | 84542 | 8456 | 8458 | 84624 | 84639 | 84668 | 84676 | 84690 | 84692 | 84699 | 84700 | 84734 | 84740 | 84796 | 84804 | 84852 | 84911 | 84929 | 84935 | 84951 | 84958 | 84983 | 85329 | 85416 | 85439 | 85474 | 8608 | 8609 | 8647 | 8708 | 8710 | 8755 | 8784 | 8807 | 8814 | 8823 | 8844 | 8909 | 8928 | 8973 | 89792 | 89866 | 89932 | 90199 | 90204 | 90233 | 90321 | 9034 | 9050 | 90736 | 9080 | 9099 | 91010 | 91056 | 9113 | 91181 | 9125 | 913 | 91319 | 91369 | 9152 | 9154 | 91754 | 9180 | 91801 | 91894 | 9205 | 9210 | 9214 | 922 | 92283 | 9230 | 926 | 9260 | 92610 | 9265 | 92691 | 92960 | 9317 | 93550 | 93624 | 9394 | 93979 | 94081 | 94104 | 94235 | 9457 | 9487 | 9509 | 9534 | 958 | 9597 | 9617 | 9650 | 9687 | 970 | 9704 | 9712 | 9807 | 9819 | 9848 | 9902 | 993 | 9942 | 9965 | 9966 | 9967 | 9980

Cluster 5:
100009676 | 100127983 | 100128327 | 100128816 | 100128893 | 100130071 | 100130264 | 100130285 | 100131131 | 100131454 | 100131496 | 100132356 | 100141515 | 10018 | 10023 | 10047 | 100506742 | 10053 | 10060 | 10110 | 10117 | 1014 | 10154 | 10158 | 1017 | 10205 | 10223 | 10232 | 10249 | 10255 | 10260 | 10324 | 10331 | 10332 | 10333 | 10343 | 10344 | 10345 | 1038 | 10389 | 10499 | 10526 | 10555 | 10578 | 10666 | 107 | 10733 | 10803 | 10846 | 10858 | 10876 | 10880 | 10886 | 10893 | 10913 | 10984 | 11012 | 11035 | 11055 | 11063 | 11064 | 11074 | 11081 | 11086 | 11101 | 11119 | 11126 | 11144 | 11185 | 11194 | 11202 | 11223 | 11249 | 11254 | 11261 | 112802 | 112840 | 112937 | 113277 | 1135 | 1136 | 1137 | 113746 | 1143 | 1144 | 114571 | 114770 | 114780 | 114781 | 114784 | 114821 | 114881 | 114905 | 114960 | 114987 | 115019 | 115350 | 115352 | 115701 | 115761 | 116039 | 116113 | 116211 | 116285 | 117153 | 117157 | 117159 | 117195 | 117286 | 117532 | 117583 | 118490 | 118856 | 1193 | 119467 | 120863 | 121273 | 121278 | 121504 | 122046 | 122481 | 122553 | 122664 | 122665 | 1235 | 1238 | 124221 | 124590 | 124637 | 124817 | 125061 | 125972 | 126147 | 126549 | 126669 | 126695 | 126859 | 127 | 127343 | 127534 | 127540 | 127731 | 128025 | 128153 | 128488 | 1286 | 128817 | 128854 | 128861 | 128876 | 129685 | 130106 | 130576 | 130749 | 130827 | 130888 | 130940 | 130951 | 131177 | 131368 | 131578 | 131669 | 132228 | 132241 | 132243 | 132851 | 132864 | 132989 | 133396 | 133688 | 133690 | 134288 | 134353 | 134957 | 135656 | 1358 | 135932 | 1361 | 136227 | 136242 | 136263 | 136306 | 136332 | 1364 | 137075 | 137797 | 137902 | 138240 | 138255 | 138724 | 139081 | 139170 | 139378 | 140432 | 140453 | 140458 | 140462 | 140564 | 140628 | 140685 | 140701 | 140803 | 140828 | 140836 | 140838 | 140894 | 1412 | 1419 | 1420 | 142679 | 142686 | 142689 | 143458 | 143630 | 143888 | 144108 | 144406 | 144577 | 145264 | 145282 | 145501 | 145645 | 145864 | 146057 | 146177 | 146434 | 146512 | 146722 | 146862 | 147323 | 147409 | 147700 | 147872 | 147948 | 147968 | 148 | 148066 | 148103 | 148345 | 148362 | 148534 | 148870 | 149018 | 1493 | 149428 | 149461 | 149647 | 149685 | 149954 | 150082 | 150135 | 150142 | 150197 | 150280 | 150465 | 150590 | 150962 | 151 | 151112 | 151477 | 151871 | 151888 | 1521 | 152195 | 152273 | 152405 | 152518 | 152573 | 152940 | 153020 | 153218 | 153579 | 153643 | 153657 | 154075 | 154214 | 154313 | 1544 | 154664 | 155006 | 155051 | 155465 | 157695 | 157724 | 157753 | 157869 | 157983 | 158158 | 158219 | 1584 | 158584 | 158809 | 158833 | 159091 | 159371 | 160140 | 160287 | 160418 | 162083 | 162461 | 162466 | 162514 | 162517 | 162540 | 162605 | 163071 | 163486 | 1636 | 163778 | 163786 | 164237 | 164656 | 165055 | 165100 | 165530 | 165631 | 165679 | 1671 | 168451 | 168620 | 168975 | 169026 | 169270 | 170063 | 170482 | 170572 | 170689 | 170691 | 170692 | 170825 | 171019 | 1757 | 1761 | 1768 | 1769 | 1773 | 1787 | 1828 | 1849 | 185 | 189 | 192111 | 1944 | 1946 | 195828 | 196074 | 196472 | 196740 | 196913 | 199221 | 199704 | 199713 | 199870 | 200132 | 200159 | 200403 | 200504 | 200523 | 200558 | 200931 | 201292 | 201501 | 201516 | 202 | 2028 | 202915 | 203100 | 203111 | 203414 | 2035 | 2038 | 2042 | 2046 | 2050 | 2069 | 207107 | 2072 | 208 | 2081 | 2099 | 2120 | 2122 | 2130 | 2150 | 2155 | 2159 | 2169 | 219623 | 219833 | 219875 | 220004 | 220965 | 221079 | 221223 | 221393 | 221409 | 221481 | 221613 | 221756 | 2218 | 221935 | 222 | 222171 | 222236 | 222659 | 222826 | 222894 | 222950 | 223117 | 2252 | 2264 | 2268 | 22797 | 22821 | 22844 | 22853 | 22869 | 229 | 22904 | 22925 | 22995 | 23101 | 23116 | 2312 | 23145 | 23200 | 2323 | 23240 | 2329 | 23354 | 23509 | 2352 | 23520 | 23553 | 23569 | 23627 | 23746 | 245711 | 245802 | 245973 | 246176 | 2491 | 252884 | 253639 | 253714 | 253738 | 253769 | 253832 | 253935 | 254187 | 254394 | 254439 | 254773 | 255082 | 255324 | 255394 | 255403 | 255631 | 256051 | 256329 | 256764 | 256933 | 256957 | 257144 | 257169 | 257240 | 257313 | 25775 | 25794 | 25818 | 259283 | 25939 | 26005 | 260429 | 26054 | 26070 | 26085 | 26148 | 26150 | 26191 | 26223 | 26231 | 26281 | 26290 | 26291 | 26298 | 2638 | 2641 | 2645 | 2646 | 26548 | 2661 | 266743 | 2672 | 2674 | 2695 | 2702 | 2703 | 27085 | 27094 | 27101 | 27115 | 27121 | 27160 | 27177 | 27190 | 27199 | 27231 | 27287 | 27330 | 2738 | 2741 | 2776 | 280664 | 2822 | 282890 | 283150 | 283208 | 283238 | 283284 | 283385 | 283521 | 283537 | 283579 | 283677 | 283687 | 283726 | 284040 | 284071 | 284161 | 284306 | 284323 | 284353 | 284355 | 284367 | 284382 | 284418 | 284440 | 2845 | 284525 | 284697 | 284805 | 284836 | 284904 | 2850 | 285154 | 285175 | 285220 | 285268 | 285598 | 285605 | 285676 | 285696 | 285740 | 2859 | 286 | 286002 | 286103 | 286128 | 286256 | 2864 | 287015 | 2899 | 2904 | 29113 | 29125 | 2917 | 2925 | 2969 | 2977 | 2984 | 29895 | 29943 | 29974 | 29999 | 3026 | 3034 | 3036 | 3038 | 3077 | 30811 | 3082 | 3083 | 30837 | 3090 | 31 | 3111 | 3117 | 314 | 3172 | 3174 | 319100 | 3195 | 32 | 3211 | 3213 | 3218 | 3242 | 3248 | 3263 | 3339 | 3340 | 335 | 3362 | 337 | 338321 | 338324 | 338440 | 338661 | 338675 | 338773 | 338809 | 339105 | 339416 | 339501 | 339834 | 339883 | 339965 | 340146 | 341359 | 342527 | 3440 | 3455 | 3456 | 345895 | 346 | 347404 | 347732 | 347735 | 347741 | 348303 | 348654 | 350 | 353144 | 353189 | 353274 | 356 | 3567 | 3592 | 360023 | 3620 | 3645 | 3676 | 3698 | 3702 | 374286 | 3743 | 374454 | 374739 | 375057 | 375189 | 3753 | 375346 | 375444 | 375449 | 3756 | 377007 | 378825 | 3792 | 3849 | 3859 | 3868 | 387718 | 387914 | 388335 | 3885 | 388507 | 388969 | 389075 | 389118 | 389161 | 389376 | 389610 | 389941 | 390928 | 391059 | 3918 | 392255 | 3929 | 3950 | 3953 | 3960 | 3982 | 399671 | 400629 | 400668 | 4007 | 400720 | 400794 | 400926 | 401052 | 401145 | 4012 | 401494 | 401551 | 402415 | 4025 | 405754 | 4065 | 4068 | 4081 | 411 | 4121 | 4146 | 415 | 4151 | 4153 | 4159 | 4179 | 4184 | 4222 | 4249 | 4295 | 4308 | 4312 | 431705 | 4359 | 4363 | 440073 | 440087 | 440184 | 440335 | 440356 | 440695 | 440712 | 440854 | 441601 | 4430 | 4481 | 4507 | 4610 | 4617 | 462 | 4632 | 4633 | 4648 | 4653 | 4773 | 479 | 4800 | 4803 | 4825 | 4838 | 4854 | 4868 | 4880 | 4883 | 4889 | 489 | 4908 | 4914 | 492307 | 493829 | 495 | 4993 | 4998 | 5009 | 5015 | 5047 | 50489 | 50508 | 5051 | 50511 | 50514 | 5053 | 50700 | 5075 | 50831 | 50834 | 50940 | 5105 | 51052 | 51145 | 51156 | 51168 | 51179 | 51206 | 51208 | 51213 | 51232 | 51281 | 51297 | 51301 | 51302 | 51314 | 51327 | 51334 | 51352 | 51360 | 51450 | 5148 | 5149 | 51633 | 5169 | 51703 | 5172 | 51725 | 51802 | 5212 | 5229 | 5239 | 5241 | 525 | 5255 | 5267 | 5286 | 5331 | 53345 | 5339 | 5345 | 53831 | 53836 | 53905 | 53944 | 54072 | 54084 | 54165 | 5429 | 54363 | 54429 | 54434 | 54456 | 54464 | 5449 | 54491 | 5452 | 54549 | 54557 | 54586 | 5459 | 54622 | 54626 | 54627 | 54661 | 54665 | 54716 | 5475 | 54797 | 54805 | 54809 | 54810 | 54852 | 54854 | 54855 | 54887 | 54898 | 54937 | 54943 | 54986 | 54989 | 55008 | 55027 | 55038 | 55057 | 551 | 55161 | 55175 | 55224 | 55231 | 55237 | 55267 | 55275 | 55283 | 55363 | 554207 | 55503 | 5551 | 55510 | 55511 | 5553 | 55554 | 55653 | 55668 | 5576 | 55784 | 55787 | 55908 | 5592 | 55972 | 56 | 56121 | 56130 | 56155 | 56156 | 56163 | 5617 | 5624 | 56246 | 56287 | 5644 | 5651 | 5655 | 5657 | 56605 | 56651 | 56704 | 56729 | 56751 | 56891 | 56918 | 56938 | 56956 | 570 | 57016 | 57057 | 57096 | 57105 | 57113 | 57127 | 57151 | 57157 | 5737 | 57455 | 57463 | 57477 | 57481 | 57491 | 57520 | 57554 | 57582 | 57604 | 57608 | 57626 | 57684 | 57697 | 57700 | 5775 | 5778 | 57829 | 5784 | 58158 | 5817 | 58497 | 58511 | 58524 | 5858 | 5874 | 5890 | 58985 | 5929 | 5933 | 59335 | 59341 | 5972 | 5989 | 6005 | 6007 | 60385 | 60401 | 60436 | 60494 | 6101 | 6134 | 6231 | 6239 | 6247 | 6261 | 6295 | 6331 | 6332 | 6338 | 6344 | 635 | 6356 | 6368 | 6370 | 6373 | 638 | 639 | 63943 | 640 | 6403 | 6405 | 64063 | 64066 | 64072 | 64102 | 64106 | 64129 | 6419 | 642273 | 642280 | 64240 | 6425 | 642987 | 643797 | 6439 | 64434 | 64577 | 6469 | 64693 | 646962 | 646982 | 647174 | 6474 | 64750 | 64799 | 64800 | 64806 | 64866 | 649179 | 64922 | 6493 | 6504 | 65061 | 6517 | 6519 | 6528 | 6536 | 6542 | 65986 | 65987 | 660 | 6665 | 6671 | 6688 | 6694 | 671 | 6725 | 6769 | 6779 | 6795 | 680 | 6801 | 6813 | 6819 | 6820 | 6822 | 6871 | 6884 | 6887 | 6898 | 6899 | 6906 | 6975 | 7021 | 7046 | 7047 | 7051 | 7053 | 7068 | 7069 | 7080 | 7092 | 7140 | 7161 | 7187 | 7223 | 725 | 7273 | 728 | 728485 | 728613 | 7287 | 7288 | 7293 | 729359 | 7294 | 729866 | 7348 | 7391 | 7401 | 7421 | 7429 | 7479 | 7484 | 7490 | 7499 | 7516 | 7535 | 7549 | 7567 | 7592 | 7597 | 7627 | 763 | 765 | 7745 | 776 | 7772 | 779 | 7809 | 78990 | 79345 | 79363 | 79370 | 79413 | 79444 | 79465 | 79574 | 79616 | 79669 | 79703 | 79725 | 79740 | 79741 | 79774 | 79782 | 79797 | 79800 | 79810 | 79815 | 79818 | 79822 | 79827 | 79832 | 7984 | 79840 | 79849 | 79862 | 799 | 79905 | 79906 | 79915 | 79924 | 79935 | 79955 | 80005 | 80018 | 80036 | 80054 | 80072 | 80099 | 80131 | 80224 | 80264 | 80306 | 80704 | 80739 | 8076 | 80833 | 81025 | 81030 | 81501 | 81550 | 81569 | 81607 | 81623 | 81624 | 81704 | 81792 | 81794 | 81932 | 8200 | 827 | 8320 | 8323 | 8332 | 8336 | 83473 | 83597 | 83697 | 83716 | 83850 | 83873 | 83881 | 83886 | 8390 | 83902 | 83903 | 8399 | 84063 | 84070 | 84071 | 84072 | 841 | 84100 | 84171 | 8418 | 84189 | 8419 | 84203 | 84226 | 84296 | 8435 | 84517 | 84539 | 84547 | 84561 | 84614 | 84616 | 84618 | 84623 | 84629 | 8463 | 84636 | 8464 | 84643 | 84657 | 84659 | 84660 | 84665 | 84671 | 84701 | 84733 | 84766 | 84807 | 84808 | 84839 | 84849 | 84889 | 84910 | 84919 | 8492 | 84924 | 84931 | 84945 | 84959 | 84985 | 8515 | 8516 | 85289 | 85320 | 85364 | 8538 | 85397 | 8542 | 85479 | 85509 | 85569 | 8581 | 8590 | 8602 | 8626 | 8633 | 8688 | 8689 | 8711 | 8740 | 8792 | 8793 | 8796 | 8797 | 8809 | 8817 | 8820 | 8828 | 8854 | 8856 | 8929 | 8976 | 89778 | 89822 | 89858 | 89876 | 899 | 8993 | 90019 | 90134 | 9032 | 90353 | 9052 | 9058 | 90586 | 9066 | 9075 | 90987 | 911 | 91120 | 912 | 91526 | 91662 | 91687 | 917 | 9175 | 9177 | 91828 | 9194 | 92 | 92345 | 92346 | 9242 | 92949 | 93010 | 93034 | 931 | 93107 | 9319 | 932 | 93273 | 93349 | 9340 | 9365 | 9376 | 9389 | 9402 | 94033 | 941 | 94122 | 94137 | 9427 | 9496 | 9506 | 9508 | 9545 | 957 | 9620 | 9630 | 96459 | 9668 | 9690 | 971 | 9731 | 9734 | 974 | 9744 | 9824 | 9837 | 9874 | 9880 | 9923 | 9949 | 9971 | 9990

Cluster 6:
10002 | 10003 | 10004 | 100124700 | 100128035 | 100128236 | 100128979 | 100129033 | 100131439 | 100131510 | 100169752 | 100190949 | 100507203 | 10114 | 10143 | 1015 | 10599 | 10642 | 10644 | 10663 | 10879 | 11131 | 11148 | 11219 | 11280 | 113220 | 114026 | 114836 | 116369 | 117144 | 117285 | 1179 | 120066 | 120400 | 121214 | 122258 | 122651 | 122945 | 1232 | 1237 | 124359 | 124626 | 124773 | 124961 | 125115 | 1259 | 126204 | 126402 | 127124 | 127247 | 127294 | 127435 | 127602 | 128209 | 128821 | 130013 | 130540 | 1310 | 131375 | 131450 | 133558 | 134860 | 135644 | 136541 | 136991 | 137362 | 137735 | 138199 | 139105 | 139420 | 139538 | 140456 | 140612 | 140690 | 140856 | 140870 | 1411 | 143662 | 144501 | 1448 | 146336 | 146849 | 146857 | 148231 | 148581 | 148823 | 149620 | 150353 | 150572 | 151254 | 151278 | 151658 | 1527 | 153745 | 154386 | 154865 | 155435 | 1576 | 157777 | 158035 | 158046 | 158062 | 1585 | 158830 | 159963 | 159989 | 160364 | 160728 | 160762 | 160777 | 160851 | 161003 | 161142 | 161253 | 1630 | 163259 | 163479 | 164312 | 164668 | 164684 | 165721 | 166824 | 166929 | 167 | 168400 | 168507 | 169355 | 170712 | 171389 | 171482 | 171484 | 1731 | 174 | 1750 | 1767 | 1791 | 1830 | 186 | 196477 | 196996 | 197003 | 1977 | 2001 | 201798 | 202865 | 203413 | 204219 | 204474 | 2056 | 2147 | 2165 | 218 | 219527 | 219557 | 219578 | 219938 | 219990 | 220047 | 220082 | 221395 | 221416 | 221711 | 2220 | 222546 | 222611 | 2250 | 2255 | 22798 | 22835 | 22917 | 23430 | 23439 | 23495 | 23566 | 23601 | 23624 | 24 | 245938 | 246 | 2525 | 253017 | 253512 | 254228 | 255193 | 255626 | 2559 | 2567 | 256815 | 256979 | 2570 | 257044 | 259 | 259239 | 259308 | 26077 | 26102 | 26138 | 26154 | 26228 | 2623 | 26251 | 26330 | 26526 | 266629 | 266977 | 267012 | 270 | 27004 | 27033 | 27063 | 27120 | 27133 | 27189 | 27284 | 27288 | 27290 | 283011 | 283416 | 283461 | 283471 | 283598 | 283897 | 284072 | 2841 | 284254 | 284359 | 284369 | 284403 | 284593 | 284801 | 285203 | 285313 | 285498 | 285548 | 285555 | 285735 | 285759 | 285834 | 285848 | 285852 | 285943 | 286514 | 286887 | 29018 | 29122 | 29124 | 2928 | 2986 | 2993 | 29953 | 3010 | 3080 | 30832 | 30848 | 3128 | 3198 | 3207 | 3221 | 3228 | 3233 | 3238 | 33 | 3355 | 338 | 338567 | 338785 | 339 | 339669 | 339760 | 339967 | 340273 | 340385 | 340596 | 340745 | 341019 | 343472 | 344905 | 3458 | 346007 | 346389 | 3467 | 348327 | 348738 | 348840 | 353345 | 354 | 3558 | 3578 | 3579 | 3596 | 3670 | 3713 | 3738 | 374470 | 3824 | 386653 | 388503 | 388531 | 388633 | 3888 | 388951 | 389102 | 389208 | 389400 | 389421 | 389658 | 3898 | 389840 | 3909 | 392307 | 3938 | 3981 | 400464 | 400765 | 400831 | 400866 | 400943 | 400954 | 4010 | 404203 | 404220 | 4109 | 4114 | 4115 | 4225 | 4250 | 4322 | 43849 | 440836 | 440956 | 441177 | 442247 | 4438 | 4544 | 4619 | 4626 | 4892 | 4909 | 493860 | 4956 | 4994 | 50604 | 50616 | 5077 | 50945 | 51224 | 51267 | 51270 | 51311 | 5132 | 51348 | 51384 | 5139 | 51533 | 51752 | 51806 | 5197 | 5268 | 5273 | 5275 | 5314 | 53347 | 53632 | 53841 | 53942 | 5406 | 54514 | 54546 | 5470 | 54714 | 54734 | 54766 | 54827 | 55301 | 55350 | 55370 | 55534 | 55553 | 5593 | 55971 | 5609 | 56154 | 56159 | 5618 | 5620 | 56203 | 56477 | 56547 | 56656 | 56914 | 57047 | 57111 | 57159 | 57379 | 5741 | 57713 | 58483 | 58538 | 5896 | 5897 | 59082 | 5957 | 6003 | 60482 | 619351 | 6258 | 6290 | 6337 | 6360 | 6367 | 63904 | 63978 | 6406 | 64078 | 64211 | 64241 | 642864 | 644076 | 64518 | 64581 | 646024 | 646976 | 647215 | 65217 | 653581 | 6538 | 6550 | 6564 | 6677 | 6692 | 6708 | 6718 | 675 | 6783 | 6910 | 6939 | 6943 | 6947 | 6954 | 7006 | 7166 | 7173 | 7225 | 7274 | 729121 | 729652 | 732 | 732275 | 7363 | 7399 | 7762 | 7783 | 79339 | 795 | 79689 | 79747 | 79755 | 79799 | 79820 | 79838 | 79853 | 79919 | 79931 | 79977 | 80117 | 80133 | 80157 | 80217 | 80311 | 80316 | 80323 | 80350 | 80380 | 8074 | 80763 | 81027 | 81494 | 81626 | 81851 | 8288 | 83394 | 83417 | 83446 | 83447 | 83550 | 83661 | 83741 | 83844 | 83851 | 83887 | 83983 | 84075 | 84109 | 84174 | 84451 | 845 | 84560 | 84654 | 84677 | 84694 | 84850 | 84953 | 8513 | 85290 | 8530 | 85302 | 85438 | 8556 | 8712 | 8735 | 8749 | 8838 | 8852 | 89777 | 89869 | 89870 | 89882 | 89885 | 8989 | 8994 | 8995 | 90070 | 9022 | 90342 | 9119 | 9173 | 9182 | 9271 | 92747 | 92935 | 9313 | 9333 | 93408 | 9350 | 93517 | 93661 | 9398 | 940 | 9407 | 9495 | 9573 | 959 | 9626 | 9629 | 9834 | 9963

Cluster 7:
10010 | 100127886 | 100128843 | 100129455 | 100129726 | 10055 | 10061 | 10066 | 10111 | 10140 | 10155 | 10197 | 10227 | 10229 | 10318 | 10365 | 10393 | 10454 | 10484 | 10513 | 10514 | 10525 | 10558 | 10559 | 10613 | 10668 | 10672 | 10691 | 10717 | 10765 | 10802 | 10919 | 10929 | 10939 | 11018 | 11040 | 11056 | 11059 | 11079 | 11093 | 11100 | 11102 | 11103 | 11108 | 11132 | 11153 | 11162 | 11193 | 112 | 11212 | 11218 | 112464 | 11257 | 11282 | 112939 | 112950 | 113174 | 113178 | 113189 | 113419 | 113444 | 113612 | 114034 | 114799 | 115024 | 115196 | 115209 | 115353 | 115704 | 115708 | 116535 | 117194 | 117584 | 119392 | 1201 | 121536 | 123720 | 123811 | 124540 | 125476 | 125965 | 126119 | 126308 | 126393 | 127253 | 128674 | 128869 | 129080 | 129787 | 130355 | 130916 | 131601 | 132158 | 132789 | 134553 | 1355 | 136051 | 136853 | 137492 | 138162 | 1398 | 140459 | 140467 | 1418 | 1445 | 144699 | 144717 | 145270 | 146923 | 147686 | 147965 | 148156 | 1482 | 148254 | 149371 | 1497 | 149986 | 151011 | 151648 | 1519 | 152579 | 153443 | 153733 | 155054 | 155382 | 1611 | 1613 | 161823 | 163081 | 163732 | 1649 | 170685 | 1723 | 1774 | 1797 | 1811 | 1845 | 1852 | 1855 | 1857 | 1874 | 1875 | 1876 | 1952 | 196446 | 196500 | 1998 | 199834 | 200010 | 2006 | 201266 | 2015 | 202018 | 203547 | 205251 | 2068 | 210 | 211 | 2188 | 22 | 2201 | 221477 | 221830 | 221937 | 222487 | 2239 | 22823 | 22839 | 2288 | 22944 | 23008 | 2301 | 23011 | 23013 | 23019 | 23071 | 23112 | 23132 | 23133 | 23151 | 23160 | 23173 | 23174 | 23190 | 23197 | 23216 | 23253 | 23283 | 23315 | 23353 | 23357 | 23386 | 23400 | 23415 | 2342 | 23468 | 23474 | 23549 | 23558 | 23583 | 23587 | 23603 | 23609 | 2516 | 257357 | 25776 | 25802 | 2581 | 25821 | 25903 | 25913 | 259287 | 25948 | 25970 | 25973 | 25978 | 25998 | 26003 | 26053 | 26057 | 26060 | 26064 | 26092 | 26099 | 26130 | 26133 | 26145 | 26147 | 26175 | 26225 | 26229 | 26263 | 26276 | 26284 | 26515 | 26523 | 26579 | 266727 | 26952 | 26953 | 26984 | 26996 | 27010 | 27034 | 27099 | 27161 | 27246 | 27252 | 27292 | 27300 | 27315 | 27343 | 2739 | 283254 | 283345 | 283450 | 283489 | 284312 | 284451 | 284459 | 284695 | 285195 | 285464 | 285590 | 285908 | 286075 | 286148 | 286336 | 2872 | 28977 | 28987 | 28988 | 2901 | 29068 | 2932 | 2960 | 2961 | 2965 | 2972 | 2976 | 29890 | 29919 | 29960 | 3004 | 3030 | 3061 | 3098 | 3177 | 3298 | 338674 | 339123 | 340547 | 342574 | 347853 | 3482 | 348235 | 3569 | 3601 | 3661 | 3716 | 373156 | 3741 | 3779 | 3801 | 382 | 387316 | 388722 | 3985 | 399818 | 402573 | 4084 | 4126 | 4149 | 414926 | 4154 | 4200 | 4234 | 4291 | 4296 | 4302 | 4324 | 4342 | 4352 | 440104 | 440275 | 441911 | 442092 | 4520 | 4548 | 4595 | 4600 | 4615 | 4636 | 4650 | 4664 | 4665 | 4799 | 4802 | 4846 | 4879 | 488 | 4976 | 49855 | 49856 | 5000 | 5001 | 5010 | 5045 | 50485 | 5054 | 5073 | 5078 | 50999 | 51001 | 51004 | 51008 | 51012 | 51018 | 5106 | 51067 | 51099 | 51102 | 51117 | 51171 | 51172 | 5119 | 5127 | 51291 | 51307 | 51320 | 51367 | 51490 | 51499 | 51566 | 51626 | 51652 | 51654 | 51665 | 51667 | 51692 | 51742 | 5193 | 5208 | 5309 | 53349 | 5361 | 53630 | 541565 | 54441 | 54505 | 54511 | 54532 | 54535 | 5458 | 54585 | 54587 | 54617 | 5467 | 54680 | 54707 | 54764 | 54778 | 5478 | 54800 | 54819 | 54847 | 54867 | 54890 | 54904 | 54921 | 54942 | 54952 | 54954 | 54955 | 54961 | 54965 | 54977 | 54980 | 54982 | 54984 | 55033 | 55080 | 55081 | 55093 | 55101 | 55109 | 55128 | 55149 | 55176 | 55181 | 55201 | 55204 | 55209 | 55255 | 55269 | 55286 | 55290 | 55317 | 55331 | 55337 | 55341 | 55361 | 5549 | 55520 | 55585 | 55602 | 55611 | 55646 | 55690 | 55705 | 55709 | 55716 | 55720 | 55722 | 55734 | 55735 | 55738 | 55751 | 55758 | 55764 | 55766 | 55773 | 55776 | 55793 | 55818 | 5583 | 55893 | 56254 | 56302 | 56339 | 56344 | 56623 | 56658 | 56927 | 56936 | 56964 | 56971 | 57002 | 57109 | 57120 | 57130 | 5715 | 57162 | 57222 | 57226 | 57456 | 57461 | 57522 | 57531 | 57532 | 57542 | 57587 | 57595 | 57599 | 57609 | 57630 | 57680 | 57695 | 57721 | 57722 | 57787 | 5818 | 582 | 58487 | 58528 | 5930 | 59349 | 5986 | 5993 | 60313 | 6059 | 60625 | 60682 | 6197 | 6293 | 6294 | 6302 | 6320 | 63877 | 63935 | 64145 | 64221 | 64223 | 642475 | 64342 | 64343 | 6452 | 64754 | 64773 | 64776 | 64777 | 64778 | 64789 | 647946 | 6480 | 64860 | 64919 | 65003 | 65059 | 65062 | 65080 | 65123 | 65260 | 65265 | 6582 | 65979 | 6603 | 6651 | 6656 | 6668 | 6729 | 6745 | 6772 | 6782 | 6786 | 6787 | 6794 | 6830 | 6907 | 6916 | 6942 | 7030 | 7073 | 7204 | 7257 | 7264 | 7297 | 729993 | 7376 | 741 | 7442 | 7508 | 755 | 7726 | 7746 | 7755 | 7771 | 7799 | 79003 | 79023 | 79035 | 79050 | 79084 | 79088 | 79139 | 79142 | 79154 | 79364 | 79414 | 79443 | 79576 | 79583 | 79621 | 79637 | 79663 | 79675 | 79691 | 79706 | 79710 | 79735 | 79736 | 79748 | 79796 | 79805 | 79828 | 79868 | 79874 | 79882 | 79932 | 7994 | 80003 | 80023 | 80146 | 80179 | 80196 | 80254 | 80273 | 80332 | 8065 | 80864 | 80895 | 8111 | 81502 | 81570 | 81608 | 81611 | 8170 | 81790 | 81855 | 81929 | 8226 | 8242 | 8243 | 8269 | 8321 | 833 | 83481 | 83642 | 83746 | 83752 | 83755 | 8379 | 83877 | 83940 | 83999 | 84079 | 84181 | 84216 | 84219 | 84255 | 84269 | 84277 | 84287 | 84295 | 84314 | 84447 | 84450 | 84464 | 8450 | 84524 | 8455 | 84641 | 84680 | 84696 | 84705 | 84717 | 84720 | 84752 | 84858 | 84859 | 8487 | 84893 | 84912 | 8493 | 84934 | 84957 | 84978 | 84981 | 84993 | 85001 | 85365 | 85476 | 8554 | 85865 | 865 | 867 | 8672 | 8675 | 8677 | 8694 | 87178 | 8732 | 8760 | 8763 | 8803 | 8869 | 8880 | 8895 | 8898 | 8906 | 8930 | 8934 | 896 | 89845 | 8986 | 89894 | 90313 | 90338 | 90410 | 90826 | 90835 | 90933 | 9101 | 9110 | 9117 | 91179 | 91272 | 91298 | 9135 | 91433 | 91445 | 91664 | 9169 | 9170 | 9183 | 91875 | 9191 | 9202 | 9203 | 9209 | 92106 | 92140 | 92162 | 9219 | 9231 | 92335 | 92342 | 92400 | 92521 | 92579 | 92595 | 9287 | 930 | 93058 | 93183 | 9329 | 9341 | 94005 | 9441 | 950 | 9527 | 9529 | 9554 | 9569 | 9570 | 9572 | 961 | 9669 | 9681 | 9682 | 9686 | 9755 | 9767 | 9772 | 9779 | 9815 | 9844 | 9852 | 9870 | 9873 | 9879 | 9887 | 9900 | 9913 | 9929 | 9931 | 9943 | 996 | 9962 | 9972

Cluster 8:
10006 | 10025 | 10128 | 10294 | 10295 | 10299 | 10329 | 10450 | 10617 | 10728 | 10809 | 10817 | 1102 | 11044 | 11186 | 1119 | 11232 | 11311 | 11325 | 116143 | 116150 | 121260 | 124402 | 129642 | 130557 | 131965 | 137964 | 142940 | 146691 | 147007 | 148223 | 148423 | 153129 | 171023 | 1795 | 1939 | 1979 | 200576 | 203286 | 221302 | 221545 | 224 | 22872 | 22906 | 22930 | 22976 | 22982 | 23029 | 23032 | 23107 | 23122 | 23139 | 23161 | 23163 | 23244 | 2332 | 23324 | 23325 | 23332 | 23369 | 23383 | 23410 | 23451 | 23598 | 23633 | 23708 | 23731 | 24138 | 24144 | 2475 | 2531 | 253558 | 253782 | 254170 | 255919 | 25766 | 25871 | 25880 | 25921 | 25926 | 25962 | 25963 | 25999 | 26048 | 26278 | 2631 | 27072 | 27230 | 27241 | 27250 | 27309 | 2764 | 2800 | 2801 | 283232 | 2842 | 29110 | 29789 | 29799 | 2992 | 30827 | 3612 | 3614 | 404636 | 4076 | 4216 | 4338 | 4358 | 440026 | 4528 | 4597 | 460 | 4690 | 4756 | 4848 | 4891 | 493 | 51020 | 51060 | 51073 | 51185 | 51547 | 51634 | 51727 | 5283 | 53371 | 54093 | 54407 | 54468 | 54496 | 54602 | 54708 | 54788 | 5481 | 54872 | 54902 | 54914 | 5498 | 55002 | 55031 | 55054 | 55166 | 55172 | 55197 | 55206 | 55239 | 55246 | 55319 | 55326 | 55506 | 55591 | 55697 | 55718 | 55727 | 55780 | 55781 | 55922 | 55968 | 56478 | 5663 | 56889 | 56900 | 56902 | 57020 | 57099 | 57102 | 57117 | 5728 | 57403 | 57419 | 5747 | 57472 | 57600 | 57693 | 5770 | 5869 | 59348 | 6018 | 60684 | 6341 | 6426 | 64318 | 64428 | 64430 | 64771 | 64801 | 66000 | 6612 | 6856 | 6904 | 7109 | 7270 | 7320 | 7347 | 7403 | 7507 | 7559 | 7566 | 7570 | 7580 | 7593 | 7750 | 79027 | 79048 | 79157 | 79161 | 79571 | 79685 | 79738 | 79752 | 79842 | 79892 | 79939 | 80006 | 8031 | 8050 | 80777 | 8089 | 8125 | 8140 | 81555 | 8220 | 8295 | 83440 | 83636 | 84148 | 84196 | 84285 | 84292 | 84549 | 84693 | 8480 | 84881 | 84896 | 8550 | 8697 | 8725 | 8766 | 8799 | 8890 | 8939 | 90637 | 91603 | 9181 | 9184 | 91949 | 9255 | 9344 | 93643 | 9371 | 9374 | 94097 | 9448 | 9520 | 9639 | 9679 | 9685 | 9774 | 9776 | 9778 | 9790 | 9793 | 9810 | 9812 | 9814 | 9828 | 9846 | 9847 | 987 | 9937

Cluster 9:
100129434 | 100130890 | 100134713 | 10068 | 100859930 | 10239 | 10322 | 1041 | 10500 | 10533 | 1059 | 10695 | 10714 | 10723 | 10748 | 10761 | 10855 | 10869 | 11133 | 11160 | 11247 | 11253 | 11264 | 11284 | 113730 | 115361 | 115509 | 115861 | 116071 | 116534 | 116844 | 118471 | 1186 | 123016 | 124935 | 125111 | 125950 | 126069 | 127385 | 127943 | 1293 | 1314 | 131540 | 132430 | 132612 | 134266 | 135 | 135114 | 136259 | 1370 | 137994 | 139562 | 139596 | 140686 | 144486 | 1446 | 1459 | 145942 | 145946 | 149473 | 150677 | 150864 | 151393 | 152404 | 153090 | 154881 | 157 | 157697 | 158131 | 1583 | 158696 | 1594 | 1605 | 165829 | 166968 | 167227 | 1798 | 1805 | 1813 | 191585 | 1948 | 196463 | 197 | 200765 | 201255 | 201514 | 202500 | 203054 | 205147 | 2063 | 2074 | 212 | 215 | 2187 | 219972 | 221421 | 221468 | 222235 | 22861 | 22901 | 23038 | 23062 | 23092 | 23102 | 23129 | 2319 | 23194 | 23218 | 23259 | 23277 | 23293 | 23417 | 23450 | 23510 | 23562 | 23616 | 23641 | 23644 | 24140 | 254263 | 254428 | 255101 | 256471 | 2569 | 25809 | 25940 | 26146 | 26164 | 26707 | 2671 | 26716 | 2676 | 26985 | 27023 | 27156 | 27158 | 27159 | 27179 | 27202 | 27236 | 27245 | 27289 | 27346 | 27433 | 283160 | 2833 | 283635 | 283871 | 284001 | 284716 | 285237 | 285367 | 285962 | 285966 | 2861 | 286122 | 286365 | 29070 | 29095 | 29119 | 29946 | 29979 | 30001 | 3042 | 3170 | 3274 | 3350 | 3354 | 3359 | 338398 | 341405 | 341567 | 3416 | 347736 | 348110 | 348793 | 3561 | 3586 | 3640 | 3659 | 3678 | 372 | 3728 | 374395 | 3814 | 3857 | 3861 | 3866 | 387748 | 388536 | 3891 | 390261 | 3975 | 399948 | 401827 | 4152 | 4209 | 4297 | 4298 | 433 | 434 | 440419 | 4437 | 4542 | 4613 | 4656 | 4772 | 4793 | 4842 | 4882 | 4999 | 5005 | 50801 | 50805 | 50939 | 5097 | 51057 | 51058 | 51066 | 51069 | 51106 | 5116 | 51192 | 51204 | 51231 | 51283 | 51319 | 51368 | 51379 | 51588 | 51738 | 5192 | 5261 | 5271 | 5335 | 5371 | 54345 | 54476 | 54522 | 54536 | 54676 | 54858 | 54862 | 54892 | 54922 | 54923 | 54925 | 54978 | 55009 | 55014 | 55017 | 55034 | 55066 | 55133 | 55151 | 55167 | 55218 | 5523 | 55244 | 55257 | 55280 | 55281 | 5546 | 55508 | 55544 | 55559 | 55561 | 55576 | 55659 | 55671 | 55706 | 55715 | 55846 | 55867 | 55902 | 55954 | 56000 | 56006 | 56131 | 56413 | 56675 | 56829 | 56896 | 57171 | 57175 | 57192 | 57198 | 57492 | 57510 | 5754 | 57573 | 57594 | 57705 | 57707 | 57710 | 578 | 5830 | 59352 | 5967 | 5968 | 6002 | 6016 | 60468 | 60490 | 60685 | 6102 | 63897 | 63929 | 64147 | 64149 | 64320 | 64321 | 64748 | 64755 | 64781 | 64788 | 64795 | 6483 | 64895 | 65010 | 6502 | 65094 | 65220 | 6523 | 653 | 654466 | 6573 | 6602 | 6707 | 6737 | 6924 | 6945 | 7003 | 7033 | 7084 | 7091 | 7189 | 7201 | 7268 | 7328 | 7343 | 735 | 7379 | 7392 | 7402 | 7464 | 7693 | 7779 | 78992 | 79057 | 79073 | 79100 | 79177 | 79446 | 79591 | 79629 | 79705 | 79760 | 79777 | 79786 | 79831 | 79922 | 79981 | 80021 | 80034 | 80124 | 80207 | 80256 | 8034 | 80351 | 80755 | 80772 | 80835 | 80856 | 81 | 81696 | 81697 | 81797 | 81850 | 81876 | 8202 | 8322 | 8328 | 83451 | 8349 | 83607 | 83737 | 83852 | 83985 | 83989 | 84140 | 84221 | 84366 | 84651 | 84838 | 84872 | 84923 | 84932 | 84961 | 85004 | 85026 | 8505 | 85285 | 8535 | 85360 | 8643 | 8660 | 8706 | 8743 | 8841 | 8861 | 8864 | 8891 | 8932 | 89887 | 90060 | 9024 | 9027 | 90381 | 90416 | 9088 | 91289 | 9149 | 9150 | 916 | 91775 | 91782 | 92070 | 9253 | 9261 | 9277 | 9290 | 942 | 9462 | 9472 | 9477 | 9562 | 9583 | 9683 | 9688 | 9695 | 9816 | 9991

Cluster 10:
10009 | 100129138 | 100192378 | 10043 | 10127 | 10139 | 10144 | 10162 | 10165 | 10196 | 10201 | 10241 | 104 | 10451 | 10464 | 10481 | 10489 | 10491 | 10495 | 10516 | 10520 | 10553 | 1062 | 10664 | 10693 | 10767 | 10785 | 10905 | 10924 | 10982 | 10999 | 11030 | 11046 | 1105 | 11052 | 11060 | 11070 | 11073 | 11082 | 11179 | 11190 | 112724 | 113 | 11329 | 113622 | 114041 | 114327 | 114826 | 114984 | 115 | 115825 | 115950 | 116068 | 116442 | 1181 | 1183 | 118424 | 120103 | 1203 | 123207 | 124045 | 124411 | 125228 | 126526 | 127933 | 1282 | 128387 | 129880 | 130507 | 130535 | 130814 | 1317 | 132160 | 132949 | 133015 | 133584 | 133619 | 134145 | 134359 | 134510 | 135154 | 135228 | 1362 | 137886 | 1384 | 1387 | 1390 | 1407 | 140831 | 143 | 1456 | 146050 | 147525 | 148523 | 148867 | 150274 | 150275 | 151194 | 151246 | 151473 | 151525 | 1520 | 152002 | 152006 | 152559 | 154743 | 156 | 157378 | 157680 | 158067 | 158427 | 159013 | 160 | 161742 | 162427 | 1633 | 164 | 165324 | 1657 | 166815 | 167153 | 167691 | 1678 | 169200 | 171425 | 171586 | 1735 | 1743 | 1839 | 1840 | 1879 | 192670 | 1968 | 1997 | 199746 | 200728 | 201163 | 201595 | 201965 | 202052 | 203522 | 2048 | 2057 | 2060 | 2065 | 2131 | 2135 | 2199 | 220064 | 220359 | 221120 | 221336 | 22955 | 22992 | 23012 | 23049 | 23076 | 23081 | 23086 | 23141 | 23199 | 23210 | 23242 | 23268 | 23299 | 23322 | 23327 | 23341 | 23347 | 23355 | 23404 | 23421 | 23431 | 23463 | 23511 | 23530 | 23536 | 23660 | 24149 | 254065 | 255252 | 256364 | 256949 | 257218 | 2582 | 25842 | 2588 | 25966 | 26009 | 26011 | 26031 | 26039 | 26091 | 26098 | 261734 | 2618 | 26258 | 26292 | 26512 | 26524 | 26586 | 267002 | 26974 | 27037 | 27128 | 27163 | 27229 | 2730 | 27332 | 282566 | 283143 | 283596 | 283987 | 284018 | 284439 | 284611 | 285172 | 285343 | 285521 | 285753 | 2874 | 2894 | 28962 | 28990 | 2935 | 2975 | 29761 | 29843 | 29896 | 29903 | 29915 | 29985 | 3052 | 30846 | 3087 | 3104 | 313 | 317762 | 3268 | 3382 | 338692 | 339287 | 339768 | 3419 | 344595 | 3512 | 355 | 3632 | 3652 | 369 | 3696 | 374986 | 3783 | 3799 | 3839 | 388284 | 3903 | 3908 | 3911 | 3964 | 400236 | 4013 | 4026 | 403341 | 4055 | 4089 | 4139 | 4205 | 4211 | 4240 | 4301 | 4343 | 441094 | 4534 | 4603 | 4705 | 473 | 4734 | 4795 | 4810 | 4893 | 4924 | 4928 | 4957 | 4967 | 49854 | 5013 | 50717 | 50807 | 50814 | 51000 | 51026 | 51028 | 51061 | 51088 | 51095 | 51115 | 51136 | 51143 | 51163 | 51201 | 51230 | 51249 | 51260 | 51275 | 51304 | 51317 | 51343 | 51385 | 51426 | 51429 | 51434 | 51455 | 51465 | 51477 | 51496 | 51534 | 51542 | 5155 | 51601 | 51605 | 51616 | 5167 | 51715 | 51776 | 51805 | 5194 | 5198 | 5199 | 5230 | 5311 | 5320 | 53616 | 5362 | 53916 | 54020 | 54108 | 54414 | 54465 | 54482 | 54497 | 54503 | 54619 | 54700 | 5471 | 54823 | 54826 | 54828 | 54843 | 54877 | 54883 | 54920 | 54928 | 54929 | 54930 | 54934 | 54956 | 54957 | 54973 | 54974 | 54976 | 54991 | 55071 | 55095 | 55129 | 55140 | 55153 | 55171 | 55183 | 55192 | 55213 | 55225 | 55230 | 55240 | 55254 | 55278 | 55284 | 55291 | 55295 | 55300 | 55322 | 55332 | 55333 | 55346 | 55357 | 55366 | 55421 | 5557 | 55572 | 55578 | 55609 | 55623 | 5565 | 55664 | 55677 | 55691 | 55719 | 55741 | 55743 | 55745 | 55754 | 55760 | 55778 | 55808 | 5581 | 55810 | 55835 | 5584 | 55892 | 55909 | 55917 | 55937 | 5594 | 56034 | 56062 | 56063 | 5608 | 56180 | 56603 | 56919 | 56946 | 56990 | 56999 | 571 | 57103 | 57150 | 57217 | 5734 | 57460 | 57480 | 57486 | 57511 | 57523 | 57561 | 57563 | 57575 | 57579 | 57590 | 57606 | 57620 | 57634 | 57669 | 57677 | 57696 | 57703 | 57715 | 57716 | 57805 | 5788 | 5793 | 5814 | 5825 | 58490 | 58492 | 58506 | 5871 | 5892 | 59345 | 594 | 5959 | 6041 | 60485 | 6049 | 6051 | 6248 | 6286 | 63035 | 632 | 636 | 63826 | 6385 | 63899 | 63979 | 64172 | 64224 | 64324 | 64327 | 64328 | 64417 | 64418 | 64421 | 64425 | 64429 | 64718 | 64757 | 64779 | 6478 | 648 | 64844 | 64854 | 6487 | 64897 | 64921 | 65008 | 65082 | 65083 | 65084 | 653121 | 654434 | 6617 | 6672 | 6715 | 6767 | 6840 | 6861 | 6865 | 6873 | 6874 | 6948 | 7024 | 7027 | 7043 | 7072 | 7090 | 7094 | 7156 | 7165 | 7226 | 730 | 7321 | 7336 | 7355 | 7456 | 7542 | 7564 | 7694 | 7727 | 7733 | 7786 | 7862 | 7871 | 79058 | 79071 | 79149 | 79158 | 79230 | 79269 | 79573 | 79582 | 79613 | 79634 | 79657 | 79659 | 79660 | 79684 | 79694 | 79707 | 79723 | 79744 | 79813 | 79830 | 79869 | 79886 | 79903 | 79913 | 79956 | 80031 | 80095 | 8013 | 80149 | 80155 | 80173 | 80185 | 8028 | 80308 | 80313 | 80321 | 8045 | 80714 | 80728 | 80764 | 80789 | 81537 | 81544 | 81573 | 81579 | 81602 | 81614 | 81788 | 81928 | 83479 | 835 | 83732 | 83849 | 83891 | 83932 | 83939 | 84186 | 84193 | 84236 | 84260 | 8427 | 84282 | 84312 | 84335 | 84342 | 84435 | 84449 | 8445 | 84456 | 84503 | 84514 | 84516 | 84520 | 84529 | 8453 | 84678 | 84709 | 84726 | 84749 | 8479 | 84826 | 84864 | 84897 | 84964 | 84984 | 84986 | 85406 | 85452 | 85459 | 85465 | 8563 | 8621 | 8631 | 868 | 8690 | 8698 | 8731 | 8736 | 8754 | 8780 | 8794 | 8795 | 8835 | 8837 | 88455 | 8884 | 8914 | 8974 | 8975 | 90075 | 9025 | 90293 | 904 | 90411 | 9060 | 90668 | 9069 | 9107 | 91151 | 91419 | 91584 | 91694 | 91860 | 9197 | 91975 | 91978 | 9200 | 92181 | 92344 | 9236 | 92399 | 92667 | 9275 | 92912 | 9321 | 9322 | 93474 | 93587 | 9373 | 93973 | 9400 | 9410 | 95681 | 9595 | 9631 | 9651 | 9701 | 9716 | 9754 | 9825 | 9830 | 9868 | 9871 | 9898 | 9908 | 9936 | 9958 | 9960 | 9969

Cluster 11:
1007 | 10361 | 10507 | 10681 | 10804 | 10891 | 11033 | 11061 | 11069 | 11076 | 11146 | 11155 | 11187 | 112609 | 1129 | 1131 | 114798 | 115294 | 116 | 116135 | 116443 | 118672 | 121793 | 123606 | 124056 | 124923 | 126433 | 1287 | 129049 | 132112 | 132332 | 138639 | 1392 | 140733 | 1408 | 145483 | 146433 | 146664 | 146760 | 151835 | 154661 | 157574 | 158135 | 1592 | 1607 | 1612 | 162239 | 163255 | 165215 | 169981 | 1740 | 1755 | 1844 | 1902 | 195814 | 200058 | 2166 | 2172 | 2185 | 219844 | 221692 | 222223 | 2259 | 22849 | 23040 | 23108 | 23125 | 23149 | 23228 | 23305 | 23316 | 23368 | 23504 | 23551 | 23604 | 23676 | 23732 | 23768 | 2523 | 253152 | 254531 | 2562 | 256536 | 25994 | 260293 | 26249 | 26256 | 26468 | 266722 | 27165 | 27319 | 2775 | 283212 | 284427 | 284656 | 2849 | 2868 | 286827 | 2911 | 2915 | 2981 | 29942 | 306 | 30813 | 324 | 340533 | 348980 | 357 | 3631 | 3739 | 3747 | 3786 | 3815 | 3834 | 388403 | 389125 | 400120 | 4008 | 404037 | 408263 | 4185 | 43 | 4356 | 441027 | 4935 | 503542 | 5067 | 51059 | 5136 | 51555 | 51559 | 51761 | 5303 | 5337 | 5413 | 54550 | 54769 | 54806 | 54836 | 5494 | 55030 | 55068 | 55088 | 55117 | 55212 | 55314 | 55521 | 55799 | 55800 | 5601 | 5625 | 5650 | 56648 | 56899 | 57084 | 57172 | 57213 | 57221 | 57526 | 57528 | 57537 | 57559 | 57576 | 57605 | 57613 | 57662 | 57718 | 57731 | 5798 | 58 | 5801 | 58157 | 58189 | 60680 | 613212 | 6262 | 63895 | 64137 | 64168 | 64208 | 64388 | 64399 | 6441 | 644662 | 64478 | 64793 | 64881 | 653140 | 6539 | 6545 | 65981 | 66035 | 6623 | 6770 | 6785 | 7010 | 7087 | 7103 | 7170 | 7432 | 7480 | 7504 | 783 | 785 | 7881 | 793 | 79608 | 79609 | 79628 | 79739 | 79762 | 79785 | 8027 | 80823 | 83482 | 8395 | 84253 | 8452 | 84812 | 84918 | 8499 | 85352 | 85442 | 8577 | 863 | 8821 | 8839 | 89795 | 89890 | 8999 | 90249 | 90673 | 9079 | 9223 | 9254 | 9314 | 9353 | 93664 | 9397 | 9481 | 9522 | 9607 | 9628 | 9677 | 9705 | 9722 | 9746 | 9750 | 9856 | 9892 | 9896

Cluster 12:
10014 | 10062 | 10130 | 10159 | 10160 | 10236 | 10240 | 10244 | 1025 | 10284 | 10289 | 10297 | 10313 | 10424 | 10428 | 10436 | 10446 | 10490 | 10519 | 10528 | 10554 | 10598 | 10600 | 10651 | 10656 | 10671 | 10682 | 10826 | 10898 | 10902 | 10923 | 10956 | 10960 | 10961 | 10972 | 10975 | 10998 | 11041 | 1108 | 1112 | 11244 | 11258 | 11266 | 112752 | 11321 | 11334 | 113402 | 114789 | 114790 | 114793 | 114879 | 114932 | 115703 | 116228 | 116254 | 116540 | 117246 | 119032 | 122953 | 123096 | 123169 | 124790 | 125988 | 126231 | 128240 | 130074 | 14 | 140771 | 1434 | 145553 | 146 | 146705 | 1495 | 149603 | 152217 | 1528 | 153527 | 157567 | 158293 | 1627 | 1665 | 171546 | 1729 | 1738 | 1796 | 1967 | 1981 | 1983 | 2073 | 2101 | 2137 | 219402 | 219541 | 220988 | 221061 | 221184 | 221443 | 221496 | 221927 | 2232 | 2248 | 2272 | 2273 | 22808 | 22856 | 22862 | 22875 | 22887 | 22931 | 22933 | 22934 | 22938 | 22980 | 23022 | 23046 | 23059 | 23126 | 23179 | 23234 | 23245 | 23286 | 23294 | 23352 | 23378 | 23394 | 23423 | 23424 | 23446 | 23517 | 23542 | 23560 | 23580 | 23582 | 23585 | 23640 | 23654 | 246184 | 256586 | 25804 | 25817 | 25844 | 25851 | 25875 | 25898 | 25901 | 25946 | 25949 | 25953 | 25989 | 26010 | 26024 | 26058 | 26061 | 26065 | 26122 | 26151 | 26297 | 2647 | 26504 | 26574 | 26589 | 26608 | 2693 | 27043 | 27069 | 27102 | 2734 | 27352 | 27440 | 2783 | 284339 | 284996 | 286451 | 28978 | 28991 | 28998 | 29081 | 29086 | 29087 | 29098 | 29117 | 29123 | 2926 | 2937 | 2941 | 2956 | 2962 | 2963 | 29844 | 29920 | 29945 | 29956 | 30 | 3035 | 3060 | 3066 | 30836 | 3181 | 3189 | 3295 | 3300 | 3309 | 339231 | 3418 | 3476 | 347689 | 347862 | 348487 | 3566 | 3707 | 375743 | 38 | 3840 | 3927 | 401152 | 402 | 4086 | 4122 | 4189 | 4236 | 4245 | 4329 | 4335 | 4351 | 4522 | 4601 | 463 | 4660 | 4677 | 4678 | 4692 | 4698 | 4719 | 4731 | 475 | 4898 | 4931 | 493753 | 4942 | 50810 | 5092 | 5095 | 51031 | 51056 | 51071 | 51109 | 51123 | 51125 | 51126 | 51133 | 51138 | 51149 | 51161 | 51164 | 51248 | 513 | 51310 | 51318 | 51322 | 51382 | 51390 | 51399 | 51430 | 51441 | 51447 | 51506 | 51567 | 51645 | 51646 | 51650 | 51696 | 51720 | 523 | 5257 | 526 | 5270 | 5372 | 53917 | 5394 | 54107 | 5412 | 5437 | 5439 | 54461 | 54462 | 54499 | 54675 | 54681 | 5479 | 54888 | 54931 | 54946 | 54948 | 5501 | 55037 | 55048 | 55063 | 55084 | 5510 | 5511 | 55120 | 5514 | 5518 | 55222 | 5529 | 55353 | 5538 | 55660 | 55684 | 55689 | 55739 | 55746 | 55756 | 55768 | 55794 | 55824 | 55830 | 55854 | 55856 | 55860 | 55914 | 5604 | 56252 | 5626 | 56262 | 5635 | 5636 | 56616 | 56674 | 56681 | 56683 | 5681 | 5684 | 5685 | 56888 | 5691 | 56928 | 56937 | 56954 | 56965 | 56987 | 5701 | 57017 | 57038 | 5705 | 5707 | 5708 | 57106 | 57153 | 5716 | 57178 | 57231 | 57406 | 57414 | 57418 | 57473 | 57533 | 57555 | 57558 | 57661 | 57701 | 57720 | 57804 | 5781 | 5796 | 5805 | 58190 | 58478 | 58533 | 5861 | 5862 | 5863 | 5876 | 5901 | 5911 | 593 | 6046 | 60496 | 6196 | 6304 | 6311 | 6369 | 6386 | 6389 | 64121 | 64132 | 6416 | 642 | 6427 | 64319 | 64397 | 64422 | 64431 | 6451 | 64710 | 64747 | 64965 | 64968 | 64981 | 65056 | 65095 | 65125 | 6513 | 6575 | 6594 | 65990 | 6618 | 6625 | 6643 | 665 | 6732 | 677 | 6810 | 6827 | 6834 | 6867 | 6897 | 6905 | 6919 | 694 | 6944 | 7108 | 7175 | 7260 | 728294 | 7398 | 7419 | 7423 | 7461 | 7485 | 7515 | 757 | 7572 | 760 | 7692 | 7707 | 7812 | 78999 | 79002 | 79039 | 79053 | 79068 | 79096 | 7913 | 79171 | 7919 | 79447 | 79581 | 79622 | 79719 | 79798 | 79852 | 79897 | 80017 | 80194 | 80212 | 80700 | 8073 | 8079 | 808 | 8087 | 81603 | 81631 | 81688 | 81888 | 81894 | 836 | 83605 | 83658 | 83707 | 83986 | 84067 | 84191 | 84259 | 84268 | 84275 | 84324 | 84365 | 84522 | 84626 | 84811 | 84814 | 84875 | 84878 | 84908 | 84926 | 84939 | 84988 | 84991 | 84992 | 85021 | 8520 | 8527 | 85457 | 8546 | 86 | 8611 | 8813 | 8834 | 8872 | 8883 | 8910 | 8911 | 89122 | 90226 | 90332 | 9038 | 90488 | 9049 | 90780 | 9130 | 9146 | 9147 | 9167 | 9218 | 92305 | 928 | 92999 | 9324 | 9328 | 9338 | 93380 | 94107 | 9415 | 94241 | 9429 | 9512 | 955 | 9556 | 9564 | 9640 | 9673 | 9693 | 9694 | 9698 | 9766 | 9806 | 9813 | 9839 | 9854 | 9861 | 9867 | 9878 | 9895 | 9897 | 9921 | 9948 | 997

Cluster 13:
100128775 | 10105 | 10106 | 10169 | 10273 | 10330 | 10455 | 10467 | 10540 | 10574 | 10670 | 10713 | 10897 | 10899 | 10969 | 10980 | 10985 | 10987 | 10992 | 11031 | 11034 | 11054 | 113246 | 11337 | 1155 | 118813 | 1211 | 126328 | 127687 | 1340 | 140823 | 142 | 1603 | 1653 | 1783 | 1808 | 1892 | 1982 | 2109 | 2217 | 22827 | 22920 | 22921 | 22979 | 23020 | 23113 | 23204 | 23265 | 23295 | 23385 | 23479 | 23513 | 2534 | 25824 | 26000 | 26521 | 266655 | 27044 | 27243 | 27336 | 2817 | 283991 | 284184 | 2885 | 28958 | 29085 | 2931 | 3032 | 3297 | 3646 | 374291 | 375 | 3837 | 3858 | 387 | 3895 | 3992 | 4077 | 4191 | 4259 | 4357 | 441951 | 4694 | 4695 | 4716 | 4722 | 4723 | 4726 | 4728 | 4779 | 483 | 5034 | 506 | 5093 | 51079 | 51084 | 51177 | 51202 | 51218 | 51247 | 51258 | 51287 | 51300 | 51329 | 51406 | 515 | 51504 | 51596 | 51604 | 51637 | 51647 | 51660 | 5202 | 5204 | 521 | 5213 | 5214 | 522 | 527 | 533 | 5436 | 54386 | 54662 | 54940 | 54998 | 55052 | 5521 | 55352 | 55700 | 55701 | 56270 | 56654 | 56851 | 56910 | 5693 | 5694 | 56998 | 57048 | 5713 | 57580 | 5879 | 5885 | 6050 | 6138 | 6146 | 6203 | 6303 | 6428 | 644 | 64855 | 64949 | 64975 | 65993 | 6746 | 6879 | 6923 | 6925 | 7009 | 708 | 7086 | 7126 | 7155 | 7203 | 7263 | 7284 | 7295 | 7385 | 7386 | 7389 | 746 | 7763 | 7818 | 7873 | 79038 | 79144 | 7936 | 79618 | 80757 | 81542 | 81576 | 8192 | 8209 | 8237 | 8241 | 8352 | 83543 | 84064 | 84286 | 84987 | 8574 | 8662 | 8724 | 8878 | 89927 | 90480 | 91252 | 92675 | 9296 | 9326 | 9343 | 9352 | 9361 | 9403 | 9528 | 9555 | 9638 | 9775 | 9791 | 9805 | 9915 | 9927 | 9993

Cluster 14:
100130557 | 100132815 | 10256 | 10326 | 10465 | 10551 | 10887 | 10940 | 10942 | 10991 | 11255 | 1128 | 11283 | 113235 | 1140 | 115286 | 115548 | 115727 | 116337 | 123355 | 126074 | 126353 | 126668 | 126868 | 127733 | 128602 | 1359 | 138882 | 140902 | 143187 | 144110 | 145957 | 146852 | 146894 | 147111 | 147694 | 147947 | 148113 | 148932 | 149095 | 15 | 150290 | 150365 | 150372 | 1506 | 1510 | 152185 | 153562 | 155 | 1581 | 159296 | 161291 | 161424 | 162681 | 162968 | 1638 | 164091 | 1669 | 169611 | 1747 | 1909 | 200312 | 2011 | 219409 | 221188 | 222699 | 2245 | 2275 | 22848 | 22899 | 2294 | 22953 | 23250 | 23308 | 23414 | 23481 | 23584 | 23670 | 23759 | 2524 | 254251 | 255057 | 259173 | 26188 | 26239 | 26259 | 26503 | 2692 | 27076 | 27136 | 27141 | 2769 | 2831 | 284083 | 284186 | 284273 | 2843 | 284358 | 285641 | 285973 | 286077 | 29104 | 3215 | 3351 | 3508 | 3680 | 3699 | 3705 | 3791 | 3821 | 387338 | 387893 | 388591 | 3886 | 3925 | 3996 | 399909 | 401250 | 4049 | 405753 | 414 | 4158 | 4221 | 4261 | 4634 | 4640 | 4776 | 4791 | 4824 | 4923 | 4992 | 5030 | 50837 | 51198 | 51277 | 51289 | 51306 | 5143 | 51733 | 51741 | 51816 | 5319 | 5325 | 54825 | 54891 | 55040 | 55065 | 55200 | 55298 | 5539 | 5566 | 5619 | 56673 | 5711 | 5729 | 57343 | 5739 | 57474 | 57530 | 57569 | 57646 | 57649 | 57835 | 587 | 5883 | 59269 | 59336 | 59344 | 5973 | 5983 | 60678 | 6357 | 64283 | 6461 | 646719 | 6492 | 6533 | 6610 | 6689 | 6693 | 6706 | 673 | 6891 | 6927 | 7015 | 7031 | 7148 | 729956 | 732253 | 7410 | 7512 | 7706 | 7782 | 7849 | 79056 | 79087 | 7975 | 79783 | 79794 | 79962 | 79998 | 80127 | 80169 | 80263 | 80341 | 80381 | 80742 | 80774 | 81031 | 81619 | 8178 | 83449 | 83896 | 84061 | 84166 | 84168 | 84432 | 84467 | 84513 | 8468 | 84876 | 84901 | 84904 | 85481 | 8564 | 8651 | 8728 | 8811 | 8851 | 8876 | 8879 | 8940 | 89883 | 9001 | 90268 | 9057 | 9095 | 91 | 91050 | 9213 | 923 | 924 | 92714 | 9294 | 93233 | 9390 | 94059 | 94233 | 949 | 9663 | 9675 | 9700 | 9753

Cluster 15:
10024 | 10046 | 10056 | 10058 | 10084 | 10113 | 10134 | 10171 | 10199 | 10206 | 10211 | 1022 | 10280 | 10283 | 10285 | 10342 | 10382 | 10432 | 10473 | 10487 | 10488 | 10492 | 10493 | 10539 | 10585 | 10618 | 10641 | 10660 | 10775 | 10778 | 10783 | 10801 | 10844 | 10935 | 10946 | 10988 | 11017 | 11020 | 11092 | 11124 | 11152 | 11157 | 11230 | 11231 | 11235 | 113000 | 11340 | 11342 | 114049 | 114908 | 114971 | 115098 | 116541 | 116985 | 117156 | 121665 | 124044 | 124446 | 124944 | 124995 | 131474 | 132001 | 138311 | 1453 | 146542 | 147184 | 1477 | 148479 | 1487 | 161 | 1662 | 170506 | 171017 | 1736 | 1793 | 1822 | 196513 | 196743 | 197370 | 201895 | 2021 | 203260 | 2067 | 2098 | 2107 | 2108 | 2132 | 2193 | 220070 | 2207 | 22926 | 23016 | 23066 | 23140 | 23184 | 23203 | 23212 | 23215 | 23303 | 23381 | 23384 | 23395 | 23438 | 23503 | 23512 | 2356 | 23658 | 251 | 25796 | 25814 | 25847 | 2593 | 25941 | 25942 | 25996 | 26015 | 26017 | 26063 | 26086 | 26090 | 26121 | 26123 | 26168 | 26286 | 26354 | 2639 | 26502 | 26520 | 266675 | 267 | 26995 | 27079 | 27112 | 27147 | 27153 | 2717 | 27247 | 27429 | 2778 | 2782 | 283989 | 284058 | 285636 | 285958 | 286257 | 28969 | 28971 | 28976 | 28989 | 29083 | 2958 | 29927 | 29959 | 29988 | 3073 | 30833 | 3093 | 318 | 3267 | 3275 | 3308 | 334 | 339201 | 339229 | 339230 | 3420 | 359948 | 3611 | 3654 | 378 | 391 | 409 | 4141 | 4350 | 439 | 440387 | 4594 | 4733 | 4735 | 4758 | 4762 | 48 | 4841 | 4913 | 50809 | 509 | 5096 | 51013 | 51019 | 51070 | 51075 | 51094 | 51096 | 51100 | 51112 | 51114 | 51121 | 51147 | 51154 | 51157 | 51214 | 51253 | 51292 | 51293 | 51295 | 51309 | 51428 | 51497 | 51520 | 51528 | 51530 | 51548 | 51582 | 51592 | 51608 | 51614 | 5162 | 51621 | 51629 | 51631 | 51642 | 51668 | 51693 | 51699 | 51726 | 51768 | 5226 | 5245 | 5289 | 5298 | 5306 | 5333 | 53343 | 5393 | 5431 | 5434 | 5442 | 5447 | 54480 | 54512 | 54531 | 54534 | 54584 | 54623 | 54663 | 54849 | 54865 | 54870 | 54881 | 54915 | 54949 | 54951 | 54971 | 54987 | 5499 | 550 | 55011 | 55090 | 55100 | 55116 | 55144 | 55177 | 55250 | 553115 | 55324 | 55325 | 55343 | 55379 | 55501 | 55573 | 55621 | 55625 | 55626 | 55644 | 55666 | 55699 | 55802 | 55823 | 55841 | 55848 | 55862 | 55863 | 55898 | 55905 | 55929 | 56257 | 56261 | 5690 | 56935 | 56945 | 56951 | 57128 | 57129 | 57161 | 5717 | 5719 | 5720 | 57410 | 57470 | 57479 | 57521 | 57541 | 57827 | 57862 | 5824 | 5828 | 5834 | 58477 | 58505 | 58517 | 5867 | 5914 | 5936 | 5987 | 6015 | 60314 | 60412 | 60491 | 60492 | 60560 | 617 | 6183 | 6388 | 63905 | 64067 | 64326 | 64374 | 644096 | 6448 | 6449 | 64763 | 64769 | 64792 | 64847 | 64976 | 64978 | 6524 | 65264 | 66005 | 6601 | 6621 | 6629 | 669 | 6712 | 6721 | 6723 | 6738 | 6780 | 682 | 6829 | 6835 | 6837 | 6880 | 6881 | 6883 | 7247 | 729967 | 7323 | 7332 | 7334 | 7443 | 7453 | 7555 | 7586 | 7709 | 7760 | 7803 | 7874 | 78994 | 79005 | 79065 | 79086 | 79094 | 79102 | 79145 | 79594 | 79595 | 79602 | 79631 | 79671 | 79693 | 79718 | 79720 | 79770 | 79863 | 79902 | 7993 | 80198 | 80222 | 80301 | 80324 | 80347 | 81037 | 8131 | 81556 | 81567 | 81572 | 821 | 8233 | 831 | 83546 | 83549 | 83637 | 8405 | 84173 | 84246 | 84273 | 84274 | 84304 | 84305 | 84313 | 84317 | 84328 | 84334 | 844 | 8443 | 8451 | 84545 | 84619 | 84650 | 84681 | 84975 | 8569 | 8703 | 8720 | 8726 | 8738 | 8881 | 8892 | 8907 | 89781 | 9016 | 9054 | 90809 | 9092 | 9126 | 9131 | 9158 | 9238 | 92591 | 92609 | 9266 | 92856 | 92979 | 9318 | 93487 | 9375 | 9391 | 9419 | 94234 | 9482 | 9538 | 9553 | 9559 | 9587 | 9605 | 966 | 9710 | 9726 | 988 | 994

Cluster 16:
10057 | 10179 | 10210 | 1024 | 10243 | 10363 | 10420 | 10497 | 10579 | 1060 | 10741 | 10746 | 10781 | 10806 | 10915 | 10927 | 10979 | 11019 | 1106 | 11083 | 11166 | 11196 | 11234 | 112487 | 112611 | 114044 | 114088 | 114792 | 114796 | 114882 | 115123 | 115426 | 116842 | 118 | 1182 | 118987 | 122773 | 123904 | 124808 | 129285 | 132321 | 137209 | 137695 | 143684 | 146198 | 149041 | 149420 | 150291 | 151613 | 154007 | 158234 | 158358 | 158880 | 161582 | 1628 | 1656 | 169714 | 1965 | 196527 | 196528 | 197131 | 199692 | 200894 | 200958 | 201973 | 2055 | 2182 | 219899 | 220213 | 221078 | 22850 | 22858 | 22868 | 22881 | 22889 | 22902 | 23001 | 23035 | 23168 | 23171 | 23243 | 23272 | 23279 | 23331 | 23350 | 23473 | 23527 | 23761 | 253260 | 25836 | 25852 | 259282 | 25929 | 25936 | 26959 | 27031 | 27148 | 27327 | 282809 | 284131 | 285267 | 28966 | 28999 | 29066 | 29072 | 29899 | 3096 | 3176 | 3299 | 331 | 339745 | 374868 | 3845 | 387921 | 403313 | 404201 | 4092 | 4212 | 4683 | 5007 | 50640 | 50804 | 51 | 51068 | 51105 | 51122 | 51134 | 51166 | 51167 | 51244 | 51265 | 51531 | 51593 | 5189 | 52 | 53373 | 5469 | 54756 | 54799 | 54850 | 54908 | 55015 | 55028 | 55119 | 55196 | 55219 | 55329 | 55449 | 55617 | 55634 | 55650 | 55669 | 55696 | 55703 | 55717 | 55729 | 55770 | 55870 | 56260 | 56647 | 56980 | 56985 | 56986 | 57037 | 57107 | 57147 | 57337 | 5745 | 57458 | 57482 | 57556 | 57577 | 57585 | 57597 | 57614 | 57631 | 57708 | 57711 | 58508 | 59343 | 5965 | 60481 | 60526 | 60592 | 619279 | 6198 | 6259 | 6310 | 63915 | 63920 | 64062 | 64398 | 64412 | 64645 | 6477 | 6505 | 65124 | 6541 | 65977 | 66036 | 6683 | 6760 | 6838 | 6901 | 6908 | 7181 | 7220 | 7275 | 7349 | 7352 | 7518 | 7541 | 7561 | 7626 | 7738 | 7741 | 7752 | 7776 | 79072 | 79109 | 79598 | 79677 | 79731 | 79898 | 79970 | 80135 | 80208 | 80255 | 80262 | 80745 | 81545 | 81617 | 81846 | 8301 | 83548 | 83943 | 84108 | 8411 | 84146 | 84162 | 84307 | 8446 | 8476 | 84837 | 84851 | 8500 | 85449 | 8555 | 8620 | 8776 | 8899 | 8945 | 9014 | 90639 | 90678 | 9100 | 91137 | 9129 | 91452 | 9221 | 92270 | 93627 | 94015 | 9425 | 9532 | 9692 | 9736 | 9760 | 9786 | 9924 | 9975 | 9989 | 9994

Cluster 17:
10063 | 10092 | 10102 | 10213 | 10269 | 10296 | 10328 | 10422 | 10425 | 10430 | 10440 | 10456 | 10477 | 10480 | 10552 | 10573 | 10589 | 10594 | 10632 | 10654 | 10907 | 10922 | 10938 | 10957 | 11137 | 11267 | 11273 | 11345 | 115939 | 116988 | 118487 | 119504 | 122509 | 122961 | 128308 | 1312 | 134492 | 138716 | 1499 | 155061 | 1642 | 1650 | 1652 | 200186 | 201254 | 2058 | 21 | 219927 | 221154 | 221908 | 2287 | 22937 | 23 | 23061 | 23380 | 2339 | 23412 | 23471 | 23484 | 23568 | 23673 | 246181 | 252839 | 254863 | 25793 | 25813 | 25874 | 25900 | 25912 | 25972 | 25983 | 26073 | 26100 | 26277 | 27183 | 27238 | 27257 | 27335 | 284129 | 285148 | 2879 | 28972 | 29079 | 29082 | 29101 | 29102 | 3028 | 3069 | 3281 | 3301 | 3396 | 353 | 3692 | 3831 | 4043 | 4048 | 4097 | 4292 | 4709 | 4712 | 4713 | 4718 | 4724 | 4729 | 4809 | 51024 | 51097 | 51116 | 51150 | 51188 | 51234 | 51324 | 51371 | 51433 | 51493 | 51503 | 51507 | 51510 | 51538 | 51639 | 51651 | 5184 | 5236 | 53339 | 537 | 53938 | 5438 | 5440 | 54529 | 54926 | 55004 | 55012 | 55049 | 55651 | 55658 | 55661 | 55831 | 55837 | 55967 | 55969 | 56672 | 5682 | 5683 | 5686 | 56922 | 56942 | 56943 | 5699 | 5702 | 5709 | 57136 | 57325 | 5763 | 57659 | 5768 | 5786 | 58525 | 5908 | 5917 | 5955 | 5976 | 60528 | 60626 | 6117 | 6118 | 6301 | 6396 | 64110 | 64219 | 64432 | 64743 | 64928 | 64983 | 65005 | 6520 | 6631 | 6634 | 7251 | 7311 | 7329 | 7384 | 7415 | 78988 | 79081 | 790955 | 79135 | 79169 | 7917 | 79174 | 79590 | 79672 | 79803 | 8106 | 8165 | 8239 | 83460 | 83930 | 83941 | 8409 | 84153 | 84271 | 84299 | 84320 | 84326 | 84892 | 8533 | 8625 | 8668 | 8669 | 88745 | 8888 | 8943 | 90324 | 90843 | 9114 | 91304 | 91582 | 9276 | 9367 | 9380 | 9416 | 9446 | 9550 | 9637 | 9665 | 9703 | 9727 | 9801 | 9817 | 9842 | 9987

Cluster 18:
10098 | 1012 | 10120 | 10235 | 10397 | 10423 | 10669 | 10966 | 1113 | 1114 | 112703 | 11322 | 11332 | 114794 | 116138 | 116984 | 118442 | 122416 | 126823 | 127262 | 131408 | 134548 | 134549 | 140606 | 143903 | 1465 | 148014 | 150726 | 152926 | 154796 | 163175 | 1742 | 1831 | 1850 | 192683 | 2016 | 203190 | 203259 | 2036 | 2114 | 2134 | 2161 | 219539 | 219738 | 221424 | 222166 | 2258 | 2261 | 2274 | 2281 | 22876 | 22891 | 22941 | 23017 | 23037 | 23096 | 23136 | 23154 | 23158 | 23207 | 23236 | 23241 | 23313 | 23362 | 23429 | 23466 | 23507 | 23531 | 23612 | 23646 | 246330 | 254778 | 2560 | 256472 | 25789 | 25864 | 25924 | 26056 | 26240 | 26580 | 26960 | 26999 | 27239 | 2744 | 2805 | 28232 | 283209 | 286097 | 2902 | 2905 | 29114 | 2913 | 2949 | 29965 | 30845 | 30850 | 311 | 322 | 3241 | 3251 | 340542 | 374378 | 3748 | 374875 | 3798 | 388341 | 389337 | 3914 | 40 | 4135 | 4137 | 4155 | 4208 | 4294 | 4345 | 4638 | 4674 | 4675 | 4815 | 492311 | 5019 | 5101 | 51062 | 51085 | 51090 | 51162 | 5121 | 51286 | 51299 | 51341 | 51466 | 5164 | 5168 | 51700 | 5176 | 51807 | 5305 | 54463 | 54715 | 54899 | 5502 | 5526 | 5528 | 55531 | 55643 | 55665 | 5600 | 5602 | 56895 | 56950 | 57026 | 57030 | 57062 | 57165 | 57180 | 57338 | 57415 | 57446 | 57512 | 57552 | 57596 | 57648 | 57689 | 57698 | 577 | 59338 | 59353 | 60370 | 6305 | 6327 | 64077 | 64112 | 64130 | 64405 | 64506 | 6456 | 6457 | 64834 | 64838 | 65055 | 65078 | 6585 | 65985 | 66004 | 6620 | 6792 | 6844 | 753 | 7532 | 768096 | 770 | 782 | 79041 | 79137 | 79623 | 7991 | 79929 | 80206 | 80758 | 818 | 825 | 8314 | 83719 | 83988 | 84258 | 84276 | 84302 | 84315 | 84440 | 84448 | 84628 | 85300 | 85315 | 8534 | 8618 | 8774 | 88 | 8927 | 8936 | 9215 | 9256 | 93185 | 9351 | 9399 | 9456 | 9467 | 9625 | 9670 | 9762 | 9764 | 9811 | 9853 | 9890 | 9912

Cluster 19:
10099 | 10109 | 10267 | 10383 | 10476 | 10567 | 10575 | 10724 | 10749 | 10971 | 10983 | 112483 | 11315 | 11341 | 1152 | 125144 | 1327 | 1329 | 1337 | 1347 | 1363 | 1436 | 149345 | 1655 | 1778 | 183 | 1938 | 2 | 2023 | 2286 | 23209 | 254910 | 25798 | 25823 | 26517 | 27089 | 280636 | 2823 | 2824 | 2896 | 29 | 29970 | 3094 | 3320 | 3490 | 3983 | 4082 | 4670 | 4691 | 4702 | 4708 | 4715 | 4725 | 4736 | 5052 | 50854 | 51160 | 517 | 5250 | 535 | 539 | 5516 | 55845 | 567 | 5689 | 5692 | 573 | 57348 | 5764 | 59286 | 6128 | 6156 | 6228 | 6415 | 6500 | 65009 | 6604 | 6653 | 7064 | 805 | 8099 | 84667 | 84886 | 8655 | 8664 | 8667 | 8682 | 8722 | 9045 | 9240 | 93621 | 93974 | 9612 | 9709 | 9741 | 9802

Cluster 20:
10150 | 10193 | 10302 | 10314 | 10411 | 10435 | 10444 | 10479 | 10505 | 10529 | 10771 | 108 | 10814 | 10815 | 10873 | 10947 | 113655 | 118663 | 118812 | 118924 | 1209 | 121512 | 126272 | 126731 | 127281 | 127703 | 130733 | 133022 | 135112 | 147 | 151556 | 158471 | 160518 | 166336 | 1741 | 1838 | 1995 | 20 | 200845 | 203197 | 2043 | 221002 | 221294 | 221895 | 22903 | 22911 | 23002 | 23229 | 23287 | 23314 | 23365 | 23543 | 254295 | 255374 | 25827 | 25837 | 25906 | 25945 | 26051 | 266812 | 26750 | 27067 | 27107 | 284252 | 284348 | 285527 | 286144 | 286262 | 2982 | 29967 | 326625 | 3423 | 377 | 3797 | 389072 | 400745 | 406 | 4133 | 4248 | 4355 | 440145 | 440556 | 4591 | 4644 | 4681 | 51046 | 51239 | 51333 | 5141 | 5153 | 51701 | 528 | 54472 | 5463 | 54704 | 55074 | 55342 | 55512 | 5567 | 55680 | 5575 | 55812 | 5595 | 55964 | 5638 | 56904 | 57168 | 57449 | 57478 | 57683 | 57704 | 57732 | 5799 | 5877 | 5905 | 6253 | 63925 | 64061 | 6440 | 64426 | 6543 | 6546 | 6558 | 6595 | 6812 | 6817 | 7163 | 7368 | 7390 | 7556 | 79016 | 79575 | 79641 | 79734 | 79837 | 80115 | 80210 | 80325 | 80723 | 810 | 81553 | 81565 | 83660 | 8382 | 8412 | 84264 | 84530 | 84899 | 84937 | 84938 | 8502 | 8506 | 8525 | 8536 | 8541 | 8604 | 8605 | 8646 | 8715 | 8745 | 8991 | 90634 | 9145 | 93611 | 9581 | 9644 | 9706 | 9738 | 9911 | 9922 | 9946 | 998

2.2.3 Distribuição unificada por cluster

Agora que os genes já foram agrupados por similaridade, visualiza-se o valor da informação mútua entre e intra os grupos com os genes selecionados. Para isso, será necessário discretizar as entradas para então obter as distribuições empirícas de cada gene.

https://search.r-project.org/CRAN/refmans/infotheo/html/multiinformation.html

breaks <- seq(floor(min(glioma)), ceiling(max(glioma)), by=1) # definindo intervalos para discretizar

mutualinfo_clust <- list()
for(clust in 1:length(genes.clusters)){
  cluster_disc_temp <- apply(glioma[genes.clusters[[clust]],], MARGIN = c(1,2),function(x) cut(x, breaks)) # discretizando cluster
  
  mutualinfo_clust <- append(mutualinfo_clust,multiinformation(t(cluster_disc_temp))) # Correlação Total (em nats - quando usado logaritmo natural para entropia)
}

for(i in 1:length(genes.clusters)){
  if(i==1) cat("--- Informação Mútua (em nats):\n")
  cat(paste0("\nCluster ",i,":\t ",format(round(mutualinfo_clust[[i]],2),width = 6, nsmall = 2)))
}
--- Informação Mútua (em nats):

Cluster 1:   319.02
Cluster 2:   430.44
Cluster 3:   521.60
Cluster 4:   363.50
Cluster 5:   700.54
Cluster 6:    28.82
Cluster 7:   711.06
Cluster 8:   245.52
Cluster 9:   369.32
Cluster 10:  616.31
Cluster 11:  281.62
Cluster 12:  527.31
Cluster 13:  185.21
Cluster 14:  173.84
Cluster 15:  418.11
Cluster 16:  266.11
Cluster 17:  185.53
Cluster 18:  342.20
Cluster 19:   87.61
Cluster 20:  203.14
ind_bound <- max(quantile(unlist(mutualinfo_clust),0.25)[[1]],1); ind_bound # métrica para clusters consistentes
 198.7392
clusters_consistentes <- which(mutualinfo_clust<ind_bound); clusters_consistentes # Clusters consistentes (com distribuição de intensidade luminosa similar)
  6 13 14 17 19
df_consistentes <- data.frame()
for(clust in clusters_consistentes){
  # definindo a mediana das intensidades luminosas para cada indivíduo da amostra considerando todos genes do cluster
  cluster_dist <- apply(glioma[genes.clusters[[clust]],], MARGIN=2, FUN= function(x) median(x)) 
  df_consistentes <- rbind(df_consistentes, data.frame(t(cluster_dist), row.names = clust))
}

# distribuições das intensidades lumninosas unificada por cluster
df_consistentes

2.2.4 Discriminação de características

Tem-se 3 características de cada indívíduo disponibilizadas conforme visto na seção de apresentação do conjunto de dados: gênero, idade e diagnóstico. Com isso, verifica-se se há distinção da distribuição de cada nível de tais características por cluster de genes consistentes obtidos anteriormente.

Contudo, inicialmente, é realizado um breve tratamento nos dados: para a análise dos níveis de gênero, desconsidera-se as entradas como “unknow”; além disso, categoriza-se a idade, e despreza-se os dados com entrada inválida como idade negativa; por fim, para a informação sobre diagnóstico, separa-se em 6 diagnósticos já apresentados, e em dois tipos - anaplástico e plástico.

Dessa forma, para cada cluster, tem-se a visualização da distribuição da luminosidade para esses diferentes grupos, e espera-se identificar clusters que discriminizam bem alguma característica.

valid_gender <- info |> filter(gender != "unknown") |> select(id = geoAccession, gender)

valid_age <- info |> filter(age > 0) |> select(id = geoAccession, age)
valid_age$cat_age <- cut(valid_age$age, breaks = c(0,7,14,20,30,40,50,60,70,80,100))

valid_diag <- info |> select(id = geoAccession, diagnostic) |> 
  mutate(diag_type = case_when(
    str_detect(diagnostic, "anaplastic") ~ "anaplastic",
    TRUE ~ "plastic"
    )
  )
for(clust in clusters_consistentes){
  meio_p1 <- sum(range(df_consistentes[as.character(clust),valid_gender$id]))/2
  med_p1 <- median(unlist(df_consistentes[as.character(clust),valid_gender$id]))
  leg_dir_p1 <- ifelse(med_p1<meio_p1, 1, 0)
  
  p1 <- data.frame(int_luz = unlist(df_consistentes[as.character(clust),valid_gender$id]), gender = valid_gender$gender) |> 
    ggplot(aes(x = int_luz, fill = factor(gender), color=factor(gender))) +
    geom_density(alpha = 0.2) + 
    scale_fill_manual(values = c(male="#1775ff",female="#ff371c"), name=NULL)+
    scale_color_manual(values = c(male="#1775ff",female="#ff371c"), name=NULL)+
    scale_y_continuous(expand = expansion(mult=c(0,.10)))+
    scale_x_continuous(expand= c(0,0))+
    theme(legend.position = c(leg_dir_p1,1), legend.justification = c(leg_dir_p1,1), legend.background = element_rect(fill = "#00000050"), plot.title = element_text(hjust=0), legend.key = element_rect(color="transparent", fill="transparent"))+
    labs(x= "Intensidade luminosa", y="Densidade", title="Gênero")
  
  meio_p2 <- sum(range(df_consistentes[as.character(clust),valid_age$id]))/2
  med_p2 <- median(unlist(df_consistentes[as.character(clust),valid_age$id]))
  leg_dir_p2 <- ifelse(med_p2<meio_p2, 1, 0)
  
  p2 <- data.frame(int_luz = unlist(df_consistentes[as.character(clust),valid_age$id]), age = valid_age$cat_age) |> 
    ggplot(aes(x = int_luz, fill = factor(age), color=factor(age))) +
    geom_density(alpha = 0.2)+
    scale_fill_discrete(name=NULL)+
    scale_color_discrete(name=NULL)+
    scale_y_continuous(expand = expansion(mult=c(0,.10)))+
    scale_x_continuous(expand= c(0,0))+
    theme(legend.position = c(leg_dir_p2,1), legend.justification = c(leg_dir_p2,1), legend.background = element_rect(fill = "#00000050"), plot.title = element_text(hjust=0), legend.key = element_rect(color="transparent", fill="transparent"))+
    labs(x= "Intensidade luminosa", y="Densidade", title="Faixa etária")
  
  meio_p34 <- sum(range(df_consistentes[as.character(clust),valid_diag$id]))/2
  med_p34 <- median(unlist(df_consistentes[as.character(clust),valid_diag$id]))
  leg_dir_p34 <- ifelse(med_p34<meio_p34, 1, 0)
  
  p3 <- data.frame(int_luz = unlist(df_consistentes[as.character(clust),valid_diag$id]), diagnostic = valid_diag$diagnostic) |> 
    ggplot(aes(x = int_luz, fill = factor(diagnostic), color=factor(diagnostic))) +
    geom_density(alpha = 0.2) + 
    guides(fill=guide_legend(title=NULL), color=guide_legend(title=NULL))+
    scale_y_continuous(expand = expansion(mult=c(0,.10)))+
    scale_x_continuous(expand= c(0,0))+
    theme(legend.position = c(leg_dir_p34,1), legend.justification = c(leg_dir_p34,1), legend.background = element_rect(fill = "#00000050"), plot.title = element_text(hjust=0), legend.key = element_rect(color="transparent", fill="transparent"))+
    labs(x= "Intensidade luminosa", y="Densidade", title="Diagnóstico")
  
  p4 <- data.frame(int_luz = unlist(df_consistentes[as.character(clust),valid_diag$id]), diagnostic = valid_diag$diag_type) |> 
    ggplot(aes(x = int_luz, fill = factor(diagnostic), color=factor(diagnostic))) +
    geom_density(alpha = 0.2) + 
    scale_fill_manual(values=c("anaplastic"="#5ed19f","plastic"="#c753f5"))+
    scale_color_manual(values=c("anaplastic"="#5ed19f","plastic"="#c753f5"))+
    guides(fill=guide_legend(title=NULL), color=guide_legend(title=NULL))+
    scale_y_continuous(expand = expansion(mult=c(0,.10)))+
    scale_x_continuous(expand= c(0,0))+
    theme(legend.position = c(leg_dir_p34,1), legend.justification = c(leg_dir_p34,1), legend.background = element_rect(fill = "#00000050"), plot.title = element_text(hjust=0), legend.key = element_rect(color="transparent", fill="transparent"))+
    labs(x= "Intensidade luminosa", y="Densidade", title="Tipo de diagnóstico")
  
  ptot1 <- p1+p2+plot_layout(ncol = 2) & theme(
  panel.border = element_blank(),
  panel.grid.major = element_blank(),
  panel.grid.minor = element_blank(),
  plot.background = element_rect(fill="#0c0c0c", color="#0c0c0c"),
  text = element_text(size=8)
  )
  
  ptot2 <- p3+p4+plot_layout(ncol = 2) & theme(
  panel.border = element_blank(),
  panel.grid.major = element_blank(),
  panel.grid.minor = element_blank(),
  plot.background = element_rect(fill="#0c0c0c", color="#0c0c0c"),
  text = element_text(size=8)
  )
  
  cat("<pre>--- Cluster ", clust, ":</pre>","\n")
  cat("\n")
  print(ptot1)
  print(ptot2)
  cat("\n\n")
}
--- Cluster  6 :

--- Cluster  13 :

--- Cluster  14 :

--- Cluster  17 :

--- Cluster  19 :

3 Conclusão

Repositório no GitHub